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ABSTRACT 


It had previovisly heen established that a Born^ron Kaman model 
reqin.iies interactions out to at least fifth neighbo'urs to obtain a 
good fit with the measured dispersion curves for gertaanim and silicon. 

In this dissertation a two neighbour Bom-Von Kaman model is applied 
to the diamond crjrstal. Out of the five force parameters that occur in 
the dispersion relations in the (111) direction, only one parameter is 
kept as free and is fitted with the experimental frequencies at the . 
zone boundaries iri the (lOO) and (111) directions. A good qualitative 
fit is obtained, wMch is substantially better than the fit in the case 
of gormaniura and silicon li-th a two neighbour interaction model. To 
get a better quantitative fit in case of diamond, however, interactians 
with more distant neighbours' should be taken into account. 

Secondly,- a technique based on -the use of Fourier analysis is 
suggested for' the computation of the frequency distiibution function of 
a ciystal. The .method has been applied to some two- and three-dimensional 
lattice models for which exact calculations have been done earlier. The 
present method is seen to yield good agreement with the exact resifLts , 

t 

It is show that, compaJted with the ilsUal methods, the Fourier series 
expansion method^ saves ‘fcimc and labour and has the added advantage of 
giving' a simple 'analyfeic bxpresslon for the distribution functicaii The 
first fifteen coefficients, of the Fourier teims are computed for a 
number of face-centred and;' body-c®ntred cubic lattices. 

The frequency distribution function of diamond is calculated 



using the present. technique. The specific heat of diamond- is evaluated 
in the temperature range 0°“300% and is seen to agree quite well ijith 
the experirental results. Finally the analysis of the dependence of 
the frequency and amplitude of localised mode due to an isotopic 
impurity is made . 
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CHiPTER I 

CRITICAL MALTS IS OP MODELS 

1. Introauction 

The stridy of' lattice, -vibrations is one of the f-undamental 
disciplines in the theoiy of solid state. Numerous ph37sical proper- 
ties of crystals, such as specific heat, electrical and thermal 
conducti-vity , scattering of x-rays, electrons and neutrons, infraied 
absorption, spin-lattice relaxation,- Mossbauer effect etc., require 
for their detailed ejcplanation a knov/lodge of lat-fcice dynamics of the 
crystal. New -powerful techniques, both theoretical and experimental, 
have made the study of lattice dynamics all the more interesting. 
Because of the coiiplex nature of the interatomic forces ^ about which 
detailed information is lacking, it is impossible to work out the 
lattice dynamics of the crystal in an exact way. It is, therefore, 
essential to resort to various approxiitate lattic<t dynamical models 
for the calculation of the properties of the solids. We outline in 
this chap-ber, the salient feat-Ures of such models. 

2 , Einstein and Debye- Models . 

Einstein^ was one of the first to study lattice dynamics from 

the standpoint of i-bs bearing on properties of solids . He postxila-bed 

that a simple crystal could be regarded as an assembly of atomic 

oscillators each ha-uing the same natural frequency i(^ . At low 

temperatures, this model predic-bs tho variation of the lattice heat 

capacity to be as exp (- 1iCL>/KT), K being the Boltzmann constant, 

3 

while the observed variation is of T typej however, at hi^ 

1 ' 



2 

• 

2 3 

temperat-wes the agreement is good. Debye and Born and vqa Karman 
placed lattice dynamics on a more firm footing and showed that a 
proper evaluation of the vibration frequencies of coupled atomic 
osciiP-ators leads to better results; 

Debye assumed that the lattice vibrates as if it wore an 
clastic isotropic continuum, A maximum frequency is chosen in 
order to maho the total nmber of modes equal to the number of 
classical degrees of freedom. On the basis of this model the 
specific heat at low temperatures ( T ^^12) is proportional to 
and reaches its classical value 31® at high temperatures (9^^ =’fetA.^K 
and N = total number of particles in the crystal). 

At first Debye theory had remarkable success in representing 
the varia.tion with temperature of the specific heat of a large number 
of substances. But latoi^lo\j- temperature colorimetric data definitely 
revealed that the vibration spectrum of the crystal can be very 
diffeient from the Debye bpectrum, ' The prime objection to this 
model is that it ignores the crystal anisotropy ahd the periodic 
structirce of tho crystal which causes it to bo a dispersive medium. 

4 

i^chs introduced the elastic anisotropy in the study of 

82 

specific heat of lithium, Collins has used this 'anisotropic 
continuum modol for a calculation of tho temporaturo vaniation of 
thermal expansion of cubic structuros.. This modol docs not,, 
hovrover, give a good roproscntaticai of the lattice spectrum,, The 
spectrum in this model extends to hi^or frequencies than that crm. occur 
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in a discrete lattice. 

3 , B om-von Karman Model 

3 

The Bom-von Karman model is based on an exact evaluation 
of the frequencies for the normal modes of the assembly of coupled 
atomic oscillators in a crystal. The approach involves cumbersome 
mathematics and was, therefore, overshadoixjd by the extremely simple 
Debye model. Bom-von Karman theory went relatively unnoticed till 
the experiments shoted with certainty that Debye’s model was insuffi- 
cient to explain the thermodynamic properties of solids. The plot of 
Debye chmacteristic temperature Q, for each value of C^, against T was 

not found to be a stroi^it line parallel to the T-axis, as is expected 

5 

from Debye’s theory. Blackman applied Bom-ven Karman theory to a 
linear chain, a txro dimensional square lattice and a simple cubic 
lattice and shoiArod that this theory gives rosilts more nearly like 
the observed specific heat and emphasized the use of the exact theory 
as apposed to the approxima,tc method given by Debye , while analysing 
low temperature specific heat da.ta. 

The general theory has been dc scribed by several authors 
, 6 7 8 

(Bom and Huang , Peierls , Loibfriod , and Maradudin, Montroll and 
9 . 

Mciss ). 'Je shall give here an outline of this theory and apply it 
to the case of the diamond crystal in the next chapter. 

The basic assumption on wh:' ch Bom-von Karman theory is based 
is that the nuclear motion in the lattice is determined by an effective 
policntial function ^ which depends only upon the nuclear coordinates. 
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This implies that the olcctrons are alwn,ys able to talce up a 

configuration appropriate to that of the displaced nuclei, and 

that the change of energy of electrons in the distorted crystal contri- 

10 

butes to the effective intemuclear potential. Chester has pointed 
out that the motion of the ole ctrons . will be osscntia''ly adiaba.tic, 
except for those fex^ xnth the energy close to the Fermi The 

exclusion principle discHous most of the non-adia,ba,tic transitions, 

Fc consider a three dimensional crystal composed of n unit 
cells each of x^rhich is a p.-'rallclepiped bounded by three non-coplanar 
vectors Og and a^. Each unit cell contains s atoms, fc 
the Gouilibri-um position vector of the pth unit cell with respect 
to a chosen origin by r(p) = P 2 _ + P2 ^ + Pj S-j; (P]_> P 2 P 3 

integers), whicli gives the position of one of its comers and by 
r(K) (' 'K = 1, 2, 3, 1... s ) the per, 'ti on vector of the K-th particle 
in the cell with axispect to origin of the cell. Then the position 
vector of the K-th particle in the pth unit cell shall be written as 


r (pK) = r{-p) + r (K) 


( 1 ) 


The total oner^ can be expanded as a T.yl or "ooyicr ir powers 
of the nuclear displacements u (pK) from the equilibrium configura- 
tion of tile static lattice and in the harmonic approxLmaticn it is 
given by 


+ ^ 


\ y y T- , „ 




p,p^ i,j 


(t)..(pK, pK) II. (pK) 


V" 


X u.. (p'^ K) , 


(?) 
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where is the potential energy in the equilibrixm configuration 


-0 


and u, (pK) (i = x, y, z) denotes the ith Cartesian component of the 


displacement u (pK) and 




(pK, pK) 


'1 

■>. 


4/ 


. (3) 


o- 


■^^i(pK) 2 ^ 

vihere subscript 0 signifies that, derivatives have been evaluated at the 
equilibrium positions. It may be noticed that the first order term 
in the esepansion of vanishes since in the equilibrium configuraticjn 
the force on any particle must vanish. The terms higher than the 
second order are responsible for thermal conductivity, thermal expan- 
sion, and the temperature variation of the elastic constants, as has 
been mentioned by Leibfriend and Ludwig . . The equation of motion 
of the (pK) ion of mass m^^, in ths harmonic approximation, is 

m^uL'(plO = - (pK, pK) u-(pK) (4) 

p" / 

The quantity pK) thus plays the role of the force 

constant and represents the force exerted in the ith direction on 
the atom at r(pK) when the atom r(pK) is displaced a unit distance 
in the jth direction. The periodicity of the lattice requires 

. . (pK, p/) to be S3mimetric in their indices and to. depend only 

_i- , 

on (p-p) rather than p and p^ separately: 


$ 
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JpK, pK) 


I {plC pK) = () . (p-p'^, KK) , 

j J,— X J 


(5) 


A solution of equation(4) for an independent normal mode of vibration 
of circular frequency UJ- and wave number q is 
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lii(pK) = IP(K) exp i(fl . r(pK) -LOt), 
where U^(pK) are the amplitudes. Equ.ation(4) leads to 


( 6 ) 


iVU,(K) = 21 21 


/ / 


^ i®) U.{K), 


( 7 ) 


where 


r" 


--./s 


/.h 


^± 3 ( 2 ^ KK) = Cp^.. (pK, pK) exp i a.(r(pK) - r(pK)) . 




( 8 ) 


Equatian(7) represents a set of 3s homogeneous linear equations for 
3s displacements of s atoms in the unit cell. For tne existence 
of a non- trivial solution of tloese equations j the necessaiy and 
sufficient condition is 


dot 


A 2 
D(a) - MCj 


= 0, 


(9) 




where D(g.) is a 3s x 3s matrix, lasually called the dynamical matrix, 

/ A 

whose elements are the coefficients D - - (KK, q) and M is a 3s x 3s 
diagonal matrix mth diagonal elements as m^ , m„, — , m . The 

J- eit. w 

characteristic equaticn(9) has 3s roots corresponding to a 

particular value of the wave nxmiber q. Throe of these roots, tending 
to zero as q— 3^0, are acoustic branches and the remaining 3s-3 


are optical branches. The values of q are restricted by the cyclic 

6 


boundary condition (Bom and von Karman) which requires that for 

a parallelepiped shaped crystal with its edges parallel to those 

of the unit cell and having lengths n- a- , n„a and n. s-_ 

X i. 2 ic 0 



7 


exp i q.(ii2_ai + 1^2^ + ^3^3) “ 

where n^^.rig and n^ are the mjnher of particles along the three edges 
of the -unit cell* Eqmtion ( 10 ) leads to 


q = 27S0\/n^ bq + hg/n^ bg + h^/n^ (11) 

where In , h and h_ are integers and taJse vrlxies fron zero through 
(0^- 1 ) , (^2“ 1) and 1) respectively, and b^, ^ and b^ define 
a reciprocal lattice by the equation 





( 12 ) 


In other words if the unit cell of the reciprocal lattice is subdivided 
into Uj^ngUj subccUs by dividing the three translations ■ into and n^ 

parts respectively, then (q / 2 7 T ) can only be one of the wave vectors 
connecting the origin to a point of the subdivision. For large n( = n. n-n_)_, 
the n distinct values of (£ / 2 7T ) nay be regarded as continuously and 
unif omly distributed in the reciprocal lattice cell. 


11 

Recently Lehnan et al. have proposed on axially symetric 
node! xiiich involves lesser nunber of force constants than the general 
Bom-ron Kaman nodol for the s.ane nei^^bour interaction. This 
nodel assunes that the potential energy of the crystal is a sm 
of bond-stretching and axially synnetric bond-b<aiding terns. Ihe 

first tem is proportional to the square of the conponent of the 

/ / 

relative displaconent along r (pK, pK) , the vector joining the equi- 
lihriun position of two atons at r(pK) and r(pK) , and the second tem 
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is proportional to the square of tlio relative displacement perpen- 
dicular to r(pKj pK) . The potential energy associated with the atom 
at (pK) is written as 

$2 ^ ^ [^CgCpp*', KK) 

X il. (pK, pi'6. (pK, pK))' 

+ Cg(pp^ KK) (r (pK, pK) x A r (pK, pK))^] , (13) 


where denotes a small displacement operator and C and C are the 

s B 

bond- stretching and bond-bending force constants between (pK) and 

// 

(pK) atoms. From (13) tho secular equation is obtained in the usual 
way and ono obtains 

"d (KIC^, a) = ( Skic' a(kk'', O) - a (KK, q)). ■ (14) 


The elements of A(KK', q) are given by the equation 


A. i . / 




a) = I “ 

+ C (pj)*, KK) £ 


j* / I 


2 

c(pp'> KI^) /"^q^^q. 


13 J 




exp iq.r (pK, pK) 


where C(pp, KK^ ) = Cg(pp'* KK) - CgCpp'i • 


3 

(15) 

(16) 


Tho subscripts of q in cquatian(l6) fefer to the Cartesian components 

11 11 

of q. This model has been applied to copper aluminium and 
12 

white tin with success* 


The Bom-von Karman theory involves a large number of force 
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constants. The symmeti^r elements of the crystal structure reduce the 
number of independent force constants necessarj for a general description 
of the force field appropriate to each class of neighbours . To reduce 
the number of force constants, one normally assumes short range interaction 
with the nearest and next-nearest neighbours. In the case of central 
force interactions with nearest and next-nearest neighbours, the 
secular equation for a cubic crystal contains only two force constants 
which can be related to the elastic constants by comparing the secrfLar ■ 
equation in the long wavelength limit xjith the corresponding equation 
of elasticity theory. Now for a cubic crystal, there are three elastic 
constants. For a mique determination of the force constants, there must 
be an additional rela,tion between the elastic constants* In this case, 
the model leads to the Cauchy relation C-^ C44I. Unfortunately Cauchy 

relation does not hold good for most of the metals t Experimental data 
on elastic constants show that the values of C ^ p in many cases are as 
high as two to three times the value of C44. 

Fuchs'^ attributed the failure of Cauchy relation to the presence 
of the electrons. According to him, certain energies such as ground state 
energy of the valence electron, Fermi, excliangc and the correlation 
energies of the conduction electrons, depend only on the atomic volume. 

The forces, arising from these energies, on the displaced ions cannot 
be regarded as due to the central interaction potential between the 
ions. They contribute to the compiussibdlity, but do not enter our 
consideration for strains for which there are no volume changes. The 
electrons should not, therefore, affect a purely transverse wave. The 
constants *^44 affected by the presence of the 
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electrons, beca^use they are involved in strains where the volme is 
not changed. 

Following the ideas of Fuchs, phenomenological force models 
have been introduced by the various authors . 

4 . Phenomenological Force Models CDnsideriug Electron Gas 
13 

de Launay has developed a sinplo model in which the conduction 
electrons a.re assumed to behave like an ordinary gas responding in- 
phase to the longitudinal component of the lattice waves but left 
unaffected by the transverse components. As the ^^aves travelling in 
a crystal, in general, are neither purely transverse nor purely longi- 
tuditial, tlicy arc separated into their longitudinal and transverse 
components and the electronic effect is assumed to operate in the 
longitudinal component only. The response modifies tho force components 
so that one sot is effective for longitudinal component and a different 
set for transverse component. For central force interactions with the 
nearest and next nearest neighbours, tho number of parameters that 
occur in do Lamay's model io f-ur, twice the number of parameters 
occuring in the model xd.thout the electron gas modification. These 
parameters arc evaluated from the throe elastic constants and the 
conditicn that the electron gas docs not modify the form of the secular 
equation for acoustic vjaves . The nexj relation that replaces the 
Cauchy relation is 


K = 
c 


> 


^12 ^44 


( 17 ) 



.ill 


"wiicrc K is the bullc modul’us ot tho olectron gas. From this relation 

0 

do Laxmay calculatod the number of free electrons per atom for copper 
and got a reasonable agreement with the experimental -value. By compa- 
ring the theoretical ahd experimental © - T curves for Cu> de Launay 
also observed that tho electron gas modification is important at lo-w 
temperatures ( 5°K) 

14 

Bhatia has carried out a similar investigation for cubic 
lattices , It is assumed that the forces on an ion in the lattice arise 
from a central interaction which is significant for nearest neighbours 
only and from certain energies (e.g. Fermi and exchange) which depend 
on atomic vol-ume only. Those energies are due to the free electrons. 

He is able to introduce three force constants in place of two for the 
correspondiag do Launcy model. This he achieved by a s-ui table choice 
of a term linear in displacements occuring in tho expression for 
potential.. If the exchange interaction term is dropped, the eloctrcnic 
contribution reduces to the one from do Launay’s model. 


In Bhatia’ s model 


OijCa) = {«'+oc"/3)(i-c,c,c,) + 


3 2 
4 K^a\ 


1 2 2 n 

1 + q74rr o^n2 


^ fiA 

D...(a) = s.s.c^ + 

Ij ^3./ g 1 j 


4 K a q.q. 
o ^1 0 


, 1 + K q^/4 n \^n^ 


v/here 



'..12 


^ 12 ) > ^ ^ 12 “ ^'^ 44 ) ^ 


K 


Cii - 


(b.c.c.) (18a) 


and 


A 


// 


I^ii(a) - 8 01 (3-C^Cg- CgCg- c^c^) + 4C)((2-C^Cj- C^Gg.) 


'Z n 

2 K a'^q'r 
0 


1 + ^JTo\o^ 


// 


'Di-Ca) = 4 0CS.S. + 




3 


1 + K^q^/47To^no^ 


whore 


/ „ // 

J ^1Z~ ^44^’ “ ■"^^*^11"' ‘"12“ ^*"44^' 


K. 


(2C^j_- C2_2“ ^^44^^ 


(f.c.c.) (18b) 


\vfhcrc n is tho number of oloctrons per cm and S. = sin (q.a), 
= cos (oj^a) and 2n is the lattice 
lar component of the wave numbor 2 _. 


C. = cos (a -a) and 2n is the lattice constant, q. is the jth rcctangu- 

1 ^ ri ^ 


Bhatia's theory has boon success fully applied by Jos hi and 
15 

Hernkar for a calculation of the specific heats of a number of 
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monovalent cubic metals. 

A slight variant of de Lhuhay's model has been used hj Loibfried 
and Brenig . They have employed an ideal gas in place of Sommorfold 
gas used by de Launay. However the results obtained are the same. 

17 

Becontly Sharma and Joshi have given a semiphenomenological 
model for lattice waves in metallic crystals. They visualise a meta- 
llic crystal as a lattice of bare ions imbedded in a •uniform cloud of 
electrons. The potential between ions is regarded as arising from (l) 
a central interaction between the ions , which is limited to the nearest 
and the next nearest neighbours and (2) from certain energies which 
are due to the compressibility of the conduction electrons and their 
interaction with the ions. Tho second effect is averaged over a 
¥ignor-Seitz sphere. They calculated the frequencies of normal modes 
of sodi-um in tho three symmetry directions and obtained a satisfactory 
agreement with 'bho results of neutron scattering experiments. This 
model contolns throe parameters which are related to tho experimentally 
measured clastic constants. However, the elastic constants usod by 
them have been derived by extrapolation methods . 

18 

Recently it has been pointed out by Lex that tho phenomenological 
model of dc Launay violates space group symmotiy requiremonts. In a 
b.c.c. lattice, tho points 'TT/a (1,0,0) and 7T/a (0,1,0), differing 
by a reciprocal lattice vector only, should have tho same dynamical 
matrix, but do Launay’s model gives an additional term. In -this 
respect tho models of Bhatia and that of Sharma and Joshi are also 
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not physical* Thcroforo^ the success of those models may be considered 
as a matter of chance.. 

Krebs has recently given a model that takes care of the symm- 
etry requirements by including Uiiikla.pp conditions. The effect of the 
electrons is considered through the screening parameter of the OoiiLomb 
interaction between the ions. The 5-ntcraction between closed ion shells 
is included by central interaction between the nearest and next nearest 
neighbours.. This model involves three- parameters which a^rc related 
to the three clastic constants. Krebs used the modified form of 

p A 

screening parameter given by Pines . The modification was dene in a 

21 

manner suggested bj?- Longer and Vorko . He calculated the dispersion 
curves for the b.CiCi lattices of Na and K and obtained good 
agreement with the neutron scattering results. Hoirevcr, there are 
discrepancies for Li, vrhich according to Krebs arc probably due to 
the failure of free electron approximation for k^, (the wave vector oh 
the Fermi surface). In any case, the agreement will depend on the 
choice of the screening parameter. Considering the simplicity of this 
model, it may prove very useful for considering the electron gas 
effect in mGt.als. 

5. Recent Theories: 

In the rigorous treatment of tlic lattice dynamical properties 
of metals, the main complexity originates from the interaction of 
conduction electrons with the lo.ttico. The problem is inherently a 
many body problem with a number of electrons interacting with the 
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lattice, A gonoral theory of the clectron-phonoh interaction has 

22 23 24 

been reviewed by Ziman , Chester and Sham and Ziman » Novel 

approexhes to the problem have also been disciiesod by Bardeen and 

25 20 . 26 27 

Pines , Pines , mgdal and Baym . 

Detailed studies of lattice vibrations in metals have been 
28 

carried out hy Taya . The restrictive conditions imposed in the 

theory make the applications valid for monovalent metals only. Ho 

considers the potential between ions as a sum of (1) the Coulomb 

■repulsion V . (2) the ion-ole ctron-ion interaction V„ j, which is to 

include the olcctron-olcctron interaction and (3) the sum 7^^ of the 

exchange repulsion and Van d-er Waals attraction. The first two are 

of long range nature, whereas the third one decreases very rapidly 

29 

with distance. The contribution V^ is computed by Ewald's method • 

To evaluate V^^ only the nearest neighbours arc taken into account. 

For Vju Toya used the Hartreo-Fock method in the presence of lattice 

30 

vibrations using a simplified exchange potential due to Slater . 

Toya has calculated the vibration frequencies for waves propogating 

along the symmotry directions of alkali halides and copper. In the 

case of sodim these results have been compared Tjuth the neutron scat- 

31 

toring res'ults of Woods ot al . They exhibit good agreement except 
for T^ branch in (110) direction-. Toya's theory takes care of the 
symmotrj' propertios also. It appears to bo the most fundamental one. 
The trouble , however, in its applica.tion is that it contains a number 
of olactronic parameters which are rather difficult to obtain in a^ ' 
precise manner, 
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32 

Woll cjid Kohn have discussed the theory of tho lattice 

dynamics of metals along similar lines. Their method is a modification 

33 

of Maicajina.'s quantum calculation and is essentially a Hartrec method, 
which considers only terms of the first order in both tho electron- 
electron and electron-lattice interactions. However, one cannot place 
much reliance on a calculation that ovcrloolcs the exchange effects 
and electron correlations. 

Efforts have also boon made to calculate tho atomic force 

constants from first principles-. An Ga,rly attempt in this field was 
34 

duo to Brenig . White has calculated the atomic force constants of 

35 

copper in a raoro fundamental way using Foymian’s theorem . Ho used 

36 

the Kostcr-Slater formalism , dcvolopod for localised perturbation 

in ci^Astals, to calcula.tc tho change in tho conduction oloctfoh 

chs.rgc density resulting from the infinitesimal displacement of an 

ion. Recently Cochran has interpreted the force constants of sodium 

44 

derived by Woods et al. by using the concept of psudo-potontial for 
the intcra.ction of conduction electrons with the periodic potential 
of tho lattice . The calcidation is based on the suggestion of 
Sham that tho clcctron-ion interaction potential may be regarded as 
moving rigidly with tho ions in the courso of the lattice vibrations 
and with certain limitations can be used to calculate the matrix 
elements for tho electron-phonon interaction, with the electrons 
represented by piano wa.vos. 

Experiment ally the most detailed information on lattice dynamics 
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has been obtained through the inelastic scattering of thermal and 
cold neutrons . The use of this technique has holpod in establishing, 

with reasonable accuracy, the dispersion relations for aluminixm'^^^'^^, 

42,43 , . 44 . 45 ^ . 46 . 47 - ^48 

copper ■ , sodium , magnesium , beryllium , iron , lead , 

■ V 49 X 50 , ^ n 51 . . , 52 , 53 fs 54 

niobium , tungsten , tantalum , nickel , molybdenum , W -brass, 

• 56 72 .67 , 7o' 

zinc > white tin diamond ^ germanium and silicon * 


6. Shell Model 


In the last six years some attention has been paid to the 

57 

shell model. This was originally suggested by Dick and Ovorhauser 

58 

and independently by Hanlon and Lawson to study the dielectric 

properties of alkali' halides. The first calculation for Naci was 
29 

done by Kollcman'n on the basis of Bom's theory of ionic crystals..-, 
Kcllorman treated the ions ns point charges and considered the Coulomb 
interaction between the point charges and a short range repulsiTO 
interaction between the first neighbours. The repulsive interaction 
prevents the lattice from collapsing. Kellorman-hs calculations fitted 
well with the experimental value of the frequency of the transverse 
optical mode of wave number zero. The calculated frequency distribution 
function also gave reasonable agreement with the specific heat, 

59 

Later KcU-cman-hs model liras applied to Nal by Woods ct al. 

The calculated dispersion relations wore found to bo quite - at 

variance with the experimental results obtained from neutron scattering 

data. Rj sides this, Kollorraan*'s model does not explain the dielectric 

6 

properties of the crystal. Bom and Huang pointed out that tho ratio 
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of t]n& static dielectric constant ^^to the high frequency dielectric 

constant ^^should bo )q _ q* ^^5 KGllGnaan'#»fe 

model 3 n.Glds tho value of this ratio as 3.6 as against the correct 

value 2.3 from tho moasurod dielectric constant. In this model, the 

ions are not polarizable and hence thcro is no possibility of electrons 

responding to tho applied field (which is the only contribution to the 

dielectric polarization at high frequencies). This model gives, 

therefore, the values of ^ as unity in place of nearly 2.9. 

oo 

To account for the dielectric properties of alk.ali halides, 

Dick and Overhausor introduced the shell model which has been genera.. 

60 

lisod by Cochran . The model assumes each ion to be replaced by a 

shell of charge Y and a core xd.th charge X, so that the total charge 

on the ion is (X+Y) * Each shell is connected to its owx core by an 

isotropic spring of constant K giving rise to the polarizability 
2 

(oY) /K of the free ion. In addition, there are short range coro- 
core, core-shell and shell-shell interactions. This has been shown 
in figure 1. The equations of motion arc 

A 2 AAA A AaAA. 

m U = (R + Z C Z) U + (T + Z C Y) W 

A/\.AvA 

0 = (T + Y C Z) U + (S + Y C Y) W, (19) 

A /v. A 

where R, T, S arc the matrices specifying the short range interactions 

/\ fN A 

and C is the matrix of Coulomb coefficients m,Z and Y are diagonal 
matrices representing tho mass, ionic charge, and the charge on tho 
shells respectively. 2 H h 3 ?c vectors) U is tho amplitude of tho 



F (12) 


19 



FIG. 1 THE GENERALISED SHELL MODEL DUE TO COCHRAN. 
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displaconcnts of the ions, liiilc Y W represents the Gloctrordc 

dipole ricncnt on the ions. The equations (19), with only nearest 

neighbour central short rougo forces, contain seven par-cjnctcrs . This 

nuribor increases to eight if the ionic charge is considered to be 

59 

unknown. Woods ot al. have used a simplified version of the model 

to calculate the dispersion relations for Nal. They considered only 

the negative iohs to bo polarizable and the short range interactions 

wore assumed to bo operating only through the shells; consequently 

decreasing the number of parrmotors to three which vqtq fixed by thb 

clastic constant high and low frequency dielectric constants,- 

The calculations in the syrrmetry directions were compared idth the 

results of inelastic neutron sca-ttcring experiments and with the 

predictions of the rigid ion model. Even a simplified shell model 

was scon to give much bettor agreement (Root moon - square dcvi.ation 

7^) with the experimental results than the point ion model . The 

largest deviation is near the zone boundary in tho (111) direction. 

Woods ct al. have also shown that tho discxi;prncy cannot bo removed, 

without worsening tho agreonont at other planes, by simple introduction 

of second neighbour interaction. Serious discrepancies occur for the 

longitudinal optical modes (devi..tion is as high as 20 ^ near the zone 

61 

boundary) , To remove such discrepcancios Cowley ct al. have studied' 
six different simplified versions of ■ the- more generalized Cochran's 
shell model ajid cad-culatcd tho dispersion relations for Nnl and KBr. 

The parameters -were evaluated by a least squares fit tb tin rxw.crxed 
dispersion curves in tho throe sjnanotry dirc-ctions, tho high and low 
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frequency dielectric constants end the clastic constants. “Ehc 
naxinuiii nunber of pancjictcrs wen reduced to ton by liniting the short 

m 

range interactions to the second neighbours end the ad hoc assunption 
that the short range into reactions arc through shells. 

In a nuts be 11 j the six node Is can be listed in the order of 
increasing coiriploxity as follow® . (Each includes tho features of tlie 
previous one). 

1) Only nega.tivc ions e.ro pole.rizablw end only centreJ. nce,rest 
neighbour short rongo forces are included. 

2) C'ontrel. second noe.rcst neighbour and noncontral nearest, neighbour 
forces arc included. 

3) Ionic charge is allowed to very. 

d) The polari2abilit3?' of tho positive ions is introduced assuning 
tlX’ shell charge to be tho sotac for both the ions . 

5) Positive ion polarizability is introducod in an al-tcmative vjay. 

6 

The cloctrice2 polerizability of the ion is fixed to its quoted value 

6) All the ten par .cBc tors arc allowed to vary. These include 

five short range interaction perenoters, two crystal polardzabdlitiGS , 
two short range polarizabilities and the ionic charge Z, 

Out of all these, only the sixth nodol gave good fit to the 
elastic and dielectric constants .and to tho dispersion curves in tho 
synr.i(jtry directions. HoiKJver, the difficulty with this model is that 
the short range polarizability of the positive ions turns out to bo 
negative for both liBr end Mai. On the basis of shell model it implies 
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that cither the shell, charges are positive or the spring constants 

61 

are negative, which is clearly imphysi cal. Cowley ot al. have 
suggested th.ot the inclusion of quadrupolo mononts nay improve the 
natters . 


The srne typo of cnlciiLations , without using the shell model 

^ 62,63 

conception, have also been carried out by Tolpygo and co-workers 

64,65 

and Hardy . They considered the dipole moments and took into 
account the fact that the dipole moment is dotomined by tho long 
range ns well as short range interactions. Cowley ot al. have 
shown tha.t Hardy's calculations arc essentially the saxio os obtained 
by then, except that Hardy neglected a torn that specifies short range 
interaction on, an ion by the neighbouring dipoles. Tho theory of 
Tolpygo and co-workers is also identical with the shell model theory, 
apart fron tho neglect of short range dipole-dipole forces in their 
theory. 

7. Lattice D^aiajaics of Covalent Crystals 

66 

The sitiia.tion in covalent cr3''Stals is still more obscure. Hsioch 

calculated the dispersion relations of gcnnani'un using a simple Bom- 

von Karman nodol with non-central first neighbour and central second 

neighbour interactions. I'Jhon those were conparod with the results of 

67 

inelastic neutresn scattering oxporinonts , tho agreement turned out 

to bo very poor, particxfLarly for the frequencies of tho transverse 

^68 

acoustic models where the disciopancy v/as as high as 10%, Homan 

demonstrated that by going to interactions upto fifth noi^bours end 

** * * . 
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beyond (c fifth neighbour nodol hes sixteen poranoters, but only 15 

occur in the dispersion relations in the synnotry directions), it 

is possible to obtain force nodols which fit the elastic constants end 

the expcrincntal frequencies within the oxporii;Tontal error. This 

69 

implies that forces of vovy long rnngo arc present. Lex suggested 
that tills long range nature of tho forces probably originates from 
the quadrupole-quadrupolc interaction that ojrisos duo to the distortion 
and compression of the valence electron distribution. 

60 

Cochran studied tho la.ttico dynamics of gcrmaniuin on tho 
basis of shell model, which accounts for the long range forces by 
including the interaction between electrical dipole moments induced 
on the atoms when they are* displaced from tlx'ir equilibrium positions. 
Ho considered the interactions mth nearest neighbours, involving 
six short range force constants, tho isotropic core-shell constant K 
and the shell charge Z. Tho polarizabilitjr was expressed in terms of 
the latter two parancters and tho dispersion rola.tions wore written 
in such a. wa.y th^'t only seven independent panancters occured, with 
polarizability as one of thorn. Cochran further reduced the number of 
parameters to five by putting an arbitrary rostriction th;it the ratio 
of tho X y force constant to x x force constant is tho same for the 
shoH-shcll interaction, bonding interaction between a shell and the 
surrounding atoms and tho interaction bc'twccn tho rigid atoms . Two 
of the remaining parameters vjoto fixed by the elastic constants and 
the third, was fixed by tho dielectric constant. The roinajining two 
paraiactcrs yore adjustod to fit the oxporinontal frequencies o.t the 
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cuntro and zone bomdarics in (100) and (111) directions,. The 
agrccnont in .general is quite good except for the longitudinal 
axoustic nodes of short wavelength, where the discrepancy is about ' 

14 ; 5 . 

70 

Recently Dolling tosed the shell nodol to fit the disporsuon 
curves of silicon and gomaniun. He found that a second neighbour 
model involving eleven ponanctors was required to obtain a quantitative 
fit. A genorrliscd second neighbour shell nodcl contains nine toon 
paranotors. Dolling reduced this number to eleven by imposing sonewhat 
arbitrary restrictions. 

The shell nodol has evidently proved \iscful in the case of 
alkali halides, gomaniun and silicon. The nain drawback is the large 
number of parniictcrs involved. For the NaCl structure, oven the 
first neighbour intora.ctions involve eleven parrjncters. In the diamond 
structure, the nearest neighbour interactions result in eight parameters 
and the inclusion of second nearest neighbour interactions loads to 
nineteen parameters. The number of paramo tors, of course, has been 
leducod in the various coses but not without arbitrary restrictions. 
These considerations have instigated the author to study the lattice 
dynamics of diamond on the basis of simple Bom-von Kaman model 
rather than the shell model. The dispersion relations are given in 
the next chapter. 



CHAPTER II 


LATTICE DINAMCS OF DIAMOM) 

1 , Mcxiel for Siaaond 

68 

It has previously been demonstrated by Herman that for 

germanium interactions with et least five nearest neighbours are 

required in order to get a good fit with the results from neutron 

scattering data. The same is true for silicon, as the dispersion 

71 

curves for silicon and germani-um are homologous . So many neigh- 
bours are needed probably because the transverse acoustic branches 

of germanium and silicon are very much flat near the zone boundary, 

72 

Tut for diaip.ond this is not the case . One might expect a simpler 
Rom-von Ih'-rman model to give a reasonable fit ’with the experimental 
curves. The shell model can also be tried, but the large number of 
paramotors prevents us from draxh-ng any conclusions. 

'■/ith this point in view, calculations with a second 

neighbour ■'■om-von Karman model •, uith tlie values of elastic 

consttmts and the ‘>£unan frequency as constraints, are presented here. 

In the following section, dispersion relations for diamond are deidved* 

73 74 68 

Though it hrs been done previously, * it is repeated here in 

order to explain the notation for later reference and to elaborate on 
the force constant matrices for the second neighbour interaction. 

Force Constant I-htricos for Diamond 

The diamond lattice (fig, 2) is composed of identical 
particles , lying on two equivalent interpenetrating face-centred cubic 
lattices. One of the lattices is displaced esne. quarter of the distance 


. , 25 . . 
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TIBIE 1 





^3 

k 

P 

Pi 

P2 

P3 

KK 

Fiist 
neighb Otars 
of 0 (k'= 1) 

a/2 

-a/2 

a/2 

-a/2 

a/fe 

a/2 

2 

2 

r 

2 
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0 

0 

0 

0 

21 

21 


-a/2 

.a/2 

-a/^ 

2 

3 

0 

«1 

0 

21 


a/2 

-a/2 

-a/2 

2 

4 

-1 


0 

21 


0 

0 

0 

1 

1 

0 

0 

0 

12 

First 
neighbours 
of 0 (k'=2) 

a 

a 

a 

0 

0 

a 

1 

1 

2 

3 

0 

0 

0 

1 

1 

0 

12 

12 


0 

a 

a 

1 

4 

1 

0 

0 

12 

along 'the cube diagonal of 

the other. 

Thus, 

there 

are 

two paidxLcles 

in the primitive 

unit 

coll,' 

one on 

each 

of the face-centred lattices. 


Rcferjed to a rectangular co-ordinate system with origin at a lattice 
point, the cell vectors of tte lattice are 

— I ~ ^ (0,1,1), — a (1,0,1), 3 ^ — a (l,3.,0), (20) 

-8 - 

vT^ere 2a is the lattice constant of diamond (= 3.56 x 10 cm.). The 

position vector of the equilibrium position of any particle in the 
lattice is given by equation (1). For diamond lattice the origin 
can bo chosen such that 

r(01) = (0,0,0); _r(02) = a (i.i-.iO 


( 21 ) 
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table 

Second neighbours of 0 
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P 
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a 

a 

1 

5 

0 

a 

-a 

1 

6 

a 

0 

a 

1 

7 
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0 

a 

1 

8 

a 

a 

0 

1 

9 

a 

-a 

0 

1 

lo 

0 

-a 

-a 

1 

IT 

0 

-a 

a 

1 

12 

-a 

0 

-a 

1 

13 

’E 

0 

-a 

1 

14 

-a 

•«a 

0 

1 

15 

-a 

a 

0 

1 

16 


2 (a) 

Second neighbours of 0 




^3 

K 

P 

ia 

-l-a 

-ia 

2 

5 

i-a 

-•g-a 


2 

6 

-la 

ia 

-ia 

2 

7 

3a 

"pa 

-i-a 

2 

& 

2 





-|a 


ia 

2 

9 

-|a 


ta 

2 

10 

ia 

3 

ga 

3 

2"" 

2 

11 

i-a 

|a 


2 

12 

2 

•j-a 

1 “ 

2 

13 

-la 

ia 


2 

14 

3 

-a 

2 

2 

ia 

2 

15 

2^ 

-ia 

ia 

2 

16 


Pi 

P3 

KiP. 


5 6 7 

11 12 13 

1 0 
0 -1 
1 


0 


0 

1 

0 


11 , 11 11 


8 ' 9 

14 15 

1 0 
0 0 
Jl 1 
11 11 


TABLE 2 (b) 


10 •' 

11 

12 

13 

14 

15 

16 

16 

5 

6 

7 

8 

9 

■ 10 

-1 ■ 

-1 

0 

0 

-1 

0 

1 

1 

0 

1 

-1 

0 

0 

-1 

0 

0 

-1 

0 

1 

-1 

0 

11 

U 

11 

11 

11 

11 

11 

22 

22 

22 

22 

22 

22’ 

22 


22 22 22 22 22 
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These atcans are labelled as 0 aM O' in the figure 2, Atons Ij^g 
on the same cubic lattice as 0 are labelled K = 1, and those lying 
on the same cubic lattice as 0 are labelled K ^,2. Each atom has 
four nearest neighbours arranged tetirahedrally at a distance a 43/2 
and t>;elve next nearest neighbours at a distance a 4 2, lying on 'the 
same face-centred .cubic lattice. 

Tables X, 2(a) and 2(b) give the coordinates, the corares- 

ponding numbers p and K and the integers (p. > P-r)? nearest 

X 2 o 

/ 

and next nearest neighbours of 0 and 0 , 

For calculating the dispersion relation, the force constants 
as defined bj’' equation (3) are reqtiired. Nov/ all of these force 
constants arc not independent. In fact, the number of independent 
force constants comes out to be considerably small if one exploi'ts 
the s 3 Tnmetry of the lattice. For example, the translatimal 
sjinmetry gives the relation (5), All the force constants are given 
by 3 X 3 force constant matrices between one nucleus , which can 
be chosen at the origin, and the other nuclei in the crystal; The 
force constants behave as covariant tensors of the second rank 
with respect to the coordinate transformations. In order to find 
the relations that exist between the elements of force constant 
matrices, ell those symmetry operations t^hich will bring the lattice 
into self-coincidence have to be applied. The nvmiber and nature 
of these operations are de'termined by the space group of the lattice., • 
They comprise rotations, reflections, glide planes and combinations 
of these. If a symmetry operation is expressed as a transformation 
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A. 

matrix T, then 

r (P K) = T r ( p k ) , (22) 

where r ( p k ) and r ( P K ) are two lattice points. The elements 

75 

of the matrix D (pk, p*k’ ) will then have the transformation law 
A 'TC 

D (P K, P«K») *= T D (pk, p»k») T , (23) 

v;here T is the tranpose of T and the change of indices is 
obtained from (22), 


The symmetry operations of the diamond lattice can be 

73 

described by six transformation matidces as given below. 


(i) A^: A three-fold axis of rotation, = Sg = x^, represented 

by the matrix 





0 1 
0 0 
1 0 



(24a) 


a/4j 


(ii) Ag : A centre of inversion at the point Xj_ = 
this is equivalent to a centre of inversion at (0,0,0) and a 
translation x^ = 3 ^ = Xj = a/2. The latter affects only the diange 
of label, and the ^-matrix is, therefore, that due to the inveimion 
at (0,0,0), i.e.. 



0 


-1 


0 


(24b) 
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(iii) A^: Three planes .of reflection : one is 3 ^ = giving 
the matrix 


/s 

"3 


0 10 
10 0 
0 0 1 


(24c) 


The other two planes 
obtained from and 


of re-flection = X 3 and x^ = can be 


Aj. 


(iv) A^: A rotation through TT /2 about an axis through the point 
parallel to the x^-axis foUox^red by a rotation through 
7T /2 about an axis through tiae same point parallel to the axis^ 
represented by the 'matrix 


A 




-1 0 
0 -1 



(24d) 


(v) jig: A rotation through 'Tj/S about an axis parallel to the 
x^-axis throu^ the point followed by a rotation through 

Tf /2 about an axis through the same point parallel to the x^-axis , 
with the corresponding matrix 



(24e) 



(vi) Ag : A rotation through TT about an axis parallel to the 
Xg-axLS through the point followed by a rotation through 

TT /2 about an axis through the same point parallel to the x^-axis, 
giving 


A 

T. 


0 -1 

0 0 


0 

1 



(24f) 


In addition, there are symmetry conditions given by equaticm (5). 


$ 


/ -a. 

(p, KK‘) may be represented by D where the superscript 


p denotes a particular choice of p, K and.K- , In order to use 
the formula (23) the change of indices p must also be specified. 
This is given by equation (22) and is concisely put in table 3 in 
the notation of the substitution groups. In this table by (a,b,c) 
etc. we mean that a matrix with a certain number in a cycle is to 
be derived from tlio matrix with the number following it in the 


A 


cycle by the operation T — -T i.e. = T T. Thus from we 
obtain 


A. aD a 


A P As A p ^ 

D = T D T, 
2 

Avl 


^ , p = 1, IS 


(25) 


Prom A and A_,' since D is invariant with respect to T and T^, 

1 X 
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TABLE 3 

Symmetry operations in the diamond lattice 

Syrmetry Change of label p 

operator 

A^ (1) (2 ,4,3) (5,7 ,9) (6,8,10) (15,11 ,13) (16 ,12 ,14) 

(1) (2 ,4,3) (5,7,9) (6,8,i0)(15,lT,13)(16 ,12 ,14) 

Ag (1,1) (2 ,2) (S ,3) (4,4) (5,5) (6 ,6) (7 ,7) (8 ,8) (9 ,9) 

(10,10) (n,ll) (12 ,12) (13 ,13) (14,14) (15,15) (16 ,16) 

(1) (2)(4,3)(7,5)(14,6)(12,8)(9)(16,10)(13,11)(15) 
(T) (2) (4,3) (7,5) (14,6) (12,8) (9) (16 ,10) (13 ,11) (15) 

A4 (3) (1,4, 2) (14,15, 5) (13,16 ,6) (10,12 ,7)(11, 8, 9) 

(3) (1,4,2) (14,15,5) (13,16 ,6) (10 ,l2 ,7) (11,8 ,9) 

Ag (2)(3,4,1)(10,5,13)(9,6,14)(11,7,16)(12,8,15) 

(2) (3,4,1) (10 ,5 ,13) (9,6,14) (11,7,16) (12,8,15) 

Ag (4)(2,3,1)(16,5,8)(15,6,7)(13,9,12), (14, 10,11) 

(4) (2,3,1) (16,5,8) (15,6 ,7) (13,9 ,12) (S,To,li.) 
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rience, vje may put 


(a P 

p 

P oi 


(26a) 


;''here m is the mass of the carbon atom. From A, we obtain, 


A. (n 

D = T D T , 
4 4 


This gives xis 


-(3 0^ p 

-|3 13 c4 


(26b) 


By reipeated application of A^ we find 


/v2 n 

n I 


/I 

A a4 /S 

T D T 
1 ^ 1 


P 'P 

0( ^(3 

■ p <X 


(26 c) 


^ a,2 ^ 

= I D T = - : 

1 1 


/'s, 

From A , since D is invariant with respect to T , 

U ^ 


r" 

c< 


P" 

-P 

OC 


_ P 

-P 

CiC_ 


(26d) 


A.9 A-^ a.9 A.^ , a,9 

%2' ^21' ^3"" ^3l“ ®32* 
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Therefore, ve may write 


a>9 

D 


V 

S'" 


y S 
X 


( 27 ) 


K*ow table (2b) shovra that for the second neighboiirs IL - K =1 or 2 , 

-L 

coad the integers (p^^p^^p^) for the atoms ^ and ( -^ + 6) differ 
just by a negative sign. This,' with the help of (5), pemits "US to 
write 


7p 

D 


D 


;f+6 


P = 


5 through 10 , 


(28) 


Thxis we obtain from (27) and (28), 

/!• l) 


/\15 

D 


V 


/X 

% 


% 

X 


ilso, yields. 


/S.8 

D 


i'’rom ye obtain 
5 


9 


As f\ A\ 

4 4 


T, D T, = - i 
m 




/xl5 CJlf /v8 W 
D = T_ D Tp = - 

5 5 


px V- 

V A 


s 

X 

-9' 

-*3 

X 


Tae comparison of (29) and (31) shows 


/ 


(29) 




f 


fX 


(30) 


(31) 




(32) 
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; 5 ro it may be mentioned that aoith had incorrectly ass-umed the 
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force constant matrices for the second neighbours to be sjnmiietric 
aiid hence chose ■ zero. This was, to the author's knowledge^ 

first mentioned ty Saint James”^^. 


A * 

By repeated application of A, we obtain the matrices D , 

,s5 a. 6 AlO 

D , D and D , The use of (25) and (28) gives us the rest of the 

aP 

matrices. Lastly the matrix D is obtained by noting that the 
potential energy’’ $ is invariant under rigid translations and 
rotations and honce 


^ii 

p',K« . 


Therefore , 


Dij (0, KK) 


> Dj, (pS HI'), 

pSK* 


where 'the prime on the summation sign indicates that the term 
with p' = 0, K = K’ is omitted. 


Table 4 depicts all the force constant matrices. 


TAELE 4 


Force constant matrices for the first two 
neighbours in diamond lattice. 


. .S7 


/n3 a3 
D = D = 


^■5 (^5 A 11 
D = D = D = 


ci (6 
-ji cX 
(^ -(^ cx 


A 4: /\4: 

, D = D = - 


o( -p -p 




rf^ 


<X 

A 


o< 


All 
: D = 


^7 a,7 a13 ^13 

D=D = D = D = - 


D = D = D 


x 

% 


X V 

y 




D (0,n) = D (0,2?) ^ 


r/^ ^ 

l-s x-s 

I^- V ^ 
V s 
'S X 


4-OC 


a6 a6 a12 ^12 

D = D = D =D = 


a8 a 8 

, D = D = D = D = 


X X ' S 

"■ a yiX V 

L % V ,^J 

yU, % - P 

X S 

')A /d 


aIO a 10 
D = D = 


n16 a 16 

D = D = - 


./^ " ~ s 

-■V /X 


c 


L? -S A 


10 0 

0 10 

+ L|(A+1X) 

loo 

0 10 

0 0 1 
u.- _ 


o 

o 

uL 


^ • 'Secixl a r Equation 

A 

The elements of the matrix D (g_) in (9) are obtained by 
a aimp:k, autetituticn of ^ (pK, U) , vhioh ia ^p^sented 1:^ 
fiP and given in table 5, and the corresponding coordinates of the 
different neighbours from tables 1 and 2. Then the stumiation in 
oqatuon(8) is carried out over the different citoms end the 
norron ponding elements of the dynamical matrix evaluated. For the 
model of diamond wo have chosen, those are 


%2(11) = D22.(S-2) - + X (1-C,g * /M2-CyC,-C^a^)J ■ , 
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A. 

D. 




^ 12 (^ 1 ) =‘S^p{?'2) = 4VSS -4i^S (C-C), 

x<- if. 2c y ^ z ^ X y' ’ 




^P-,(11) s. 


23^-^ - -23V-/ = - i 4 9 S^ (C^- C^) , 


/S. 


”13(11) = ”l3(»') = 4 13 S^s^ + i 4 S Sy (C^-C^) , 


A 


A 




(X 1^1 + exp - i jr|- (q^+ q^) + exp - i tl {q^+ q^) 


+ exp - i n (q^^ q^) 


%2(12) = D2q(l2) = l+cxp-ij[ (q^+ q^) 


exp - i JT (qy+ q^) - exp i /f (q^+ q^) 


023(12) = D3j,(12) = - |3 [l - exp - i ^ (q^+ q^) 

+ exp - i 7[ (qy+ q^) - exp - i n qj . 


mcl 


"0^3(12) ="03^(12) 


-PpL- exp - iff (q^+ q^) 
exp- i 7X (qy+ q^) + exp - i 7T(q2+^q^)J 


( 35 ) 
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where = cos q ), s = sxn ( 71 q ) and (Ilq7c;)is the 
x-component of the wave number q. The remaining elements are 
obtained by using the Hermitian property of the dynamical matrix. 

^ " Dispersion Curves in the Symmetry Directions 

The djmamical matrix for diamond is a 6 x 6 Hermitian matrix. 
In ^pneral, it is not possible to factorise it and one has to calcu- 
late numerically the six eigenvalues which are the noimal mode 
frequencies in the six branches . In the symmetry directions the 
dynamical matrix factorises into smaller order matrices and the 
roots of the secular equation can be mritten explicitly. Recently 

experimental dispersion curves in the (100) and (111) directions 
72 

have been reported , In both these directions the two transverse 
optic branches as i^rell as the two trrnsvorse acoustic brandies are 
degenerate. The solutions are 


(a) (100) direction: 


M Co 


M (O 


LO 


JJi } 


= 4 OC + 8 jlX ( 1 - cos TT q) 
+ 4 OC cos^ / 


(36a) 


M 00 


TO 


2 

M U) TA 


? = 4 oC + 4 ( + ^ ) ( 1 - COS TT q) 


+ 4 coe^(W/2^ + p sin(TT^/:^ . , (36b) 
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(b) (111) direction; 

2 


M CaJ 


M OJ' 


LO 


la j 


40 C+ 4 ( 2 yU.+ X + 2 D ) sin ^TT q 

? S 2 ^ 

±2 OC''(l + 3cos "Ha) + ^ ^ “ 0( )sin 7T a 




M U) TO 
2 


M U> 


TA 


(36c) 


= 4iX+4(2yU+X-l^)sin -fta 


+ 2 


^2 2 
oC (1 + 3 cos tt q)+ p ( ^ + 2oc)sin -flq 


■. 1 V 2 


(36d) 


r . E vcluation of tlie Force Constants 


For this purpose, w have four equations relating the five 

force: constants,o(. , P ’ X > and ■!) , to Raraan frequency 

end the throe elastic const, ants G. , C.,, and C , , . Ttese relations 

ll^ 12 44 

have been derived, by the long wave-length limit technique, by 


Smith 


1 0 


and 


"11 


2 a 


” 2a 


{oL + gy) , 

2 

(oC - + 4 X + 4 ) , 


= -~ (2p -C<.-4X -4/A +gD) 


( 37 ) 


2a 



• » 


The frequency shift of the first order Hanan line is 
obtained by the optical frequency at the origin. This gives 

CO = /V. 8 OC /m 


It may be noted that the force constant ^ does not 

« 

appear in the dispersion relations in the symmetrj^ directions 

and equations (37) and (38). From (37) and (38) ^ ^ 

and are expressed in terms of p and substituted in 

equations i.3C). The value of is chosen to fit the experimental 
72 

curves at the zone boundaries in the (100) and (111) directions. 

The elastic constants of dicmiond have been reported separately by 

76 • ■ 77 

Dhvagavantam and Bhimasenachar , Prince and Wooster and McSkimin 
78 

.'cnd Bond. The measurements of Prince and Wooster are comparatively 
imprecise, because the results are deduced from X-r£jy measurements 
i..a against the velocity of sound measurement by the other authors. 

The three sets of values differ qiiLte a bit from each other, as 
ciin be seen from table 5, McSkimin and Bond have discussed in 

TjfflM 5 

2 

Bias tic Constants from three different experiments , in dynes/cm 



4.4 


5.76 



detail the different sources of error in the experiment. They 
have mentioned that the error can be as high as 11,^ in the case 
of C^p due to its small value. The picture regarding the elastic 
constants appears to be quite obscure. For this reason all the 
three sets of elastic constants are tried. The value of tlie 
Raman frequency is 1332 cm“^ It has been found that 

Bhaguvontam and Bhinasenachar' s data fits best with the neutron 
scattering results ;ath the value of as ,044 x 10 dynos/cm. 

The results of this fit are shown in figure 3, The agreement^ 
in generixL, is reasonably good. The deviation 'for the LO, LA 
i,nd TO branches in both the symmetry directions is ^h.thin the 
oxper;imental error xirhich is 2-3^, For the TA branches the deviation 
is inaxirauni at ~ (100) direction and 18^ 

in the (111) direction) and decreases to less than 1% as 1. 

At this stage it tneiy be noted that the two neighbour Bom-Von 

68 

llsivacjn model for germanium results in a discrGpan.cy of about 
66^1 at the boundary of the- TA branch in the (100) direction. No 
t:hoicc 3 of (3 improves the fit. Sven a variation in the value 
of Raman frequency does not help. 

With the elastic data of Bhagavantam and Bhiraasenachar and 

6 

\-2xi vuluc of ^ as .044 x 10 dynes/cm, equations (37) and (38) 

IccKl to (in units of 10® dynos/cm) 

*. t ^ ^ *4 

CC = .157, .044/ A = ->0206, ^ = .0226, V = .0270 (39) 

It might be interesting to compare the experimental value 
of the bulk modulus K with the value obtained by our choice of the 






force constants. To calculate r we iisc the rclo-tion: 


3/K - / SaoC j (40) 

... 74 

which was suggested by Saint Jtancs on the basis of Laval's 

tlicory of elasticity. Using tte force constants from (39) and 

the elastic constants of Bhagavcntaia and Bhimasenachar, vq obtain 

from (40) 

-6 

K = 0,164 X 10 per megabar, 

wliici'. agrees vei^' well with Adam’s experimental resiilt, K = 0.16x10 

-6 

p.'r mogabnr, Hoi^ever, a different value K = O.IB x 10 per 

8l 

nogabar Vdis observod by Uillitunson . 


74 

Siiint James chooses the value of p as aero for the 

i'ol'J.ow’ng reason. If one uses the mean of the two values of 

0 

prince and Uoostcr’a data lead to a negative value of ^ /oc • 
Tills is unphysical since (y. is unambiguously dotcmiined by 
equation (38) and is positive. Instead, if one uses Adam’s value 

O 

I’or ^ positive and is very small. The value of 

P con, thcro-fore, be taken to be zero. 


Saint Jarac^s a.rgumGnt does not seem to be convincing 
because if one uses, instead of Prince and 'fooster's data, the 
otiit-r two data from table 5 ^ |3 comes out to be sufficiently 
■J,ar|.Te (/^ OC /2). The calculated dispersion curves with the 
diffcKmt sets of force constants as calculated by Saint James 
arc- found to bo in disagreoraont \n.th the experimental results. 
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The closest agreement obtained is with Bha^avantfan and Bhimasenachar's 
data md is shown in fi.gure 3. 

>•' 1 

Determination of Vi * 

1 w ’ 

Having obtained the five parametors, vjq proceed to calculate 

. It, certainly, cannot be obtained by fitting the dispersion 

relations in the symmetry directions (100) and (111), as it does 

not occur there at; all. does occur in the dispersion relations 

for (lip) direction for which no experimental da.ta exists. The 

only thing one can do is to restrict the nature of the forces. 

The gcnorilised forces can be thought of as made up of two body 

forces, three body forces, four body forces and so on. The pair 

foi'CGS arc central forces, the triplet forces are tho angular forces 

and tlic quadruplet forces are the torsional forces. By definition, 

the anfTular forces tend to resist deformation of ttao angle formed 

by cacli pair of concurrent bonds and the torsional forces tend to 

resist defomation of the dihedral angle between each pair of planes 

deh^nnined by the atoms A-B~C £ind B-G~D respectively, the atoms 

A.-P,„C~D hedng bonded in sequence, ’ifo shall assume that the second 

neighbour forces are made up of central forces , angular forces and 

torsional forces orly. The other kinds of forces for the second 

neipifiouro oie talccn to be negligible. One can show, by deriving 

13 

tiio secular equation in a ma.nnor explained by do Launey 
-nd ften comparing the secular equation for the general force 

model, that 
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D® 

central 


a9 

D 

an,<Ti.ilar 


D 

torsional 



€ 

€ 

i 

O 
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€ 

0 


0 

0 

0 


r / 

6 

/ 

€ 


_ 

€' 

e' 

2 6^ 




-4 



e' 

0 

I! 



0 


_ 0 

0 

0 


(41) 


'Hicsc individual matrices for' the central, angialar and torsional 

68 

forces have also been esdiibited by Herman . 

Therefore, ■we obtain for the restricted second neighbour 

force. 0 , 


/\9 _ /n9 ^ a9 ^ >^9 

^gf'noral " ^central angiilar torsional (42) 


From ta.ble 4 f.ind equations (41 ) and (42) ve obtain 

/ / 

%= 2 ^ and , 


so that 


6 = - X ^ 


(43) 
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Thcrofore, once we calculate X » ?) is automatically found. 

For the sake of comparison, the values of the force constants 

74 

obtained hero along x^ith those calculated by Saint J.ames aiid 
73 

Snitn arc depicted in table 6. 

TABLE 6 



Calcilations of Saint James 

Values 


;’OrC''') 

.■one tarts 

Elastic Constanto used arc dotendned 

by: 

calculated 
by Smith 
with 

Present 

calcu- 

in iini'ts 
of^ 

lO^^ipncs/ 

c:::i 

Prince and 
Wooster 

I 

MeSkimin and 
Bond 

II 

Bhagavantam & 
Ehimasenachar 

III 

constants 
of Bhagavaa- 
tam 

: lations 

IT 

oC 

0.157 

0.157 

0.157 

0.157 

0.157 

P 

0.0 

0.0 

0.0 

0.i04 

0.044 

X 

-0.0294 

-0.0163 

-0.0236 

0.0 

-0.0206 

jx 

0.0293 

0.0282 

0.0226 

0,0226 

0.0226 

D 

0.0343 

0.0312 

0.0365 

0.0226 

6.0270 

s 

0.0122 

O'. 0066 

0.0048 

0.0 

0.0103 


It nfjy be monticnod that Smith has considered a similar 
iiodt. ' wjc' opt that she restricts the second neighbour forces to only 
central forces, whereas the present model includes the .angular and 
torsi onfJ. forces ns well. For making the force central she has taken 
yU. -- ')j .For this purpose she as certs upon the 
validity of Cauchy relation established by Boro who has 
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»hown that ±n the case of central forces a cubic crystal with 

centred faces obej^s Cauchy's relation' . That relation is valid, only 

i. or a Lionatomic crystal, where every atom is a centre of ssnumetiy, 

Tho supposition of Smith, thus, amounts to neglecting the first 

neighbours T,7bich is unphysical. Though she has used the force 
o 

constant matrix, given by (41), the ar,guments used by 

l.oi' £.r'.‘ not correct. 

Tr , 83 

v'cry recently Warren has studied a different model for 

diamond, A negative charge, — e Z, situated mid ’.ray between the 

neighbouring carbon nuclei is taken to represent the charge 

dcnsi'‘iy in tho bonding direction. Tho positive charge on the carbon 

ruc:ltun rnd the .associated spherically symmetric portion of the 

electron distribution is represented by a positiw point charge eZ 

nd a point mass 11, Coulomb interactions betvjeen tte charges and 

'■herb range interactions extending up to second neighbours are 
2 

considorcid, Z is t;,jken as a parameter,. Experimental dispersion 

2 

curves in the (100) direction lead to a negative value of Z' , vfcich 
i.", clearly unphysicjfL, A least square fit with the TA(lll) branch, 
wit}), fic conr.tr.'dnt ■fch„'’.t Z is real, results in Z = 0. From this 
’•■arreo concludes that tho charge in the bonds cannot be used to 
cxplni.n the long nonge nature of forces in covalent crystals. He 
has Edso ciilculated the values of the five parameters, occuring 


in the dispersion relations in (111) direction for a second 
neighbour Bom-Von Karraan model, from the measured frequencies. 
Tho V'd-ues of the parameters so obtained are 


cxl 


« 


.1S7, 


.0236, ■)) = .0255 ^44) 


p = .0457, X = -.0222, 



These values, as can be seen from table 6, are not much different 
from the present calculations xjhere only one force constant 
is kept as a parameter. 

frequency distribution functions of diamond xdth the four 
sets, I, II, III and IV, of force constants (table 6) will be 
described in chapter ly. Use has been made of a new method, for 
t'he crlculation of frequency distribution functions, developed 
:in the next chapter. 



CHAPTER III 


CALCULATION OF FP£']UENCY SPECTRA 


1 . Int roducti on 

I'fe define the freqijencj'' distribution function, g( OO )j 

of tl'ic nonarO. modes of vibrations of crystals such that g( )d (jj 

is the fraction of frequencies in the interval- ( , CjO + CJ ) . 

2 2 

Sinilrxl.v G( OJ ) d (CU is defined to be the fraction of squared 

2 2 2 ' 

frequencies in the interval {UJ , UJ + d cO ) . The two 
functions are related by the equation, 

g( OJ ) = 2 UJ G( UJ^). (45) 

The distribution function is of fundamental importance for tlic 
unio re tending of many properties of ci^'^stals (e.g. thermodynamic 
proptortics , infrared absorption). It con be experimentally obtained 
f-or aibio crystals by using incoherent inelastic scattering of 
3 3.0V/ neutrons . Unfortunately the number of substances having 
npprc'ciublc incoherent cross~sections is limited and also the 
energy resolution of theso experiments is still too poor to reveal 
■'he f.Lnor dot/iils of g( CO ). Hi, V and Ti are the only crystals 
for vAiich the mcaEJured froq-uency spectra have been reported so 
far^’*'. In any case the dispersion curves of the normal inodes of 
vibr-etion c/in be obtained from the data on coherent scattering 
of slow neutrons and g(CO ) can, in principle, be calculated from 
these curves. In order to calculate g( CO is essential to 


. . 50 . . 
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assiaiae a force model and then fit a set of force constants to the 
dispersion curves end the elastic constants. These force constants 
give in return the value of the dynamical matrix at every point q 
in the first Brillouin zone. Usually the dispersion curves are 
measured in high symmetry directions and then analysed by the Born- 
von Karmon theory to obtain the force constants. One might question 
how adequate the physical description of the dispersicn curves is, 
as given by this nodol. This is a problem by itself and is 
iriulcvant to the present chapter which is concerned vdth a new 
re thod of computing g( tO ) . 

The methods for calculating g( Oj ) ^.re described at 

9 

length by Maradudin, Montroll and Weiss , A short review of the 
various techniques is given here so that the advantages and 
disfidv ant egos of tlic new method, that has been developed here, 
cmi bo scon, 

, F’xact LValuation 

One commonly rtsod representation of the distribution 

, 2 ^ . 86 ■ 

function G( CO ) i-s 

G ( (^^) X .Ji- S ( U)^ - U)j(k)), ('^6) 

3 s N 

j 

wterc M iB the number of unit colls anci s is the number of 
atOTS in the tinit coll* The index j is for different branches. 

The various representations of the ^ *- function (Dirac ^ function) 


388a8 
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lead, to a number of equivalent fomulae relating G( (j J^) to the 
dynamical matrix. As an application of formula (46)^ the case 
of linear lattices vath long range interactions between •(iie. atoms 
can bo considered herci 

The dispersion relation for a linear chain^ with the 
vibrations in the diioction of the chain, is 

LiJ "(k) = — i — T"” (1 - cos hk), 

PI n 

n=l 


wlioro oc ^ is the force constant for the nth neighbour. Three 

2 

diffcrcnt forms of (V are assumed: K exp( - nQ^ ), K/n 

■ n 
4 

rind K/n - Wo choose tho following representation for the delta 
function; 


^ (x) = - ^ 


substitution of (47) in (46) leads to 




(47) 


G( CO ) = — -L- iLu/v^ ^ — -T — ;; TT; 

MTT 0)^4- IC ^ 




l< 


U)^4-i€-co (k) 


Ik 




( 46 ) 


The integral in equation (4P) can oasily be calculated for the 

87 

adiclri we ha¥0 chosen. G(X) for the three cases is 
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G(X) 


(1 *. X)^ (1 - Z scch oi. /2) 


j for cy! = Ke 


-nOC 


n 


(1 - X)' 


. for oi. = — ^ 
* n 2 


L 


x^‘ (1 - }? ) 


S' 


n 


j for CX„= 


(49) 


#v ^ 

whore X = U) / uJj__ is the maxiirrum value of the 

froi^vinncy end C is the nonnalisation constant, so that 

G (X) dX = 1. (50) 


.Tt is evident from (49) that tho well known square root singularity 
in the distribution function in a linear chain persists in the 
niodeln considered, hero. 


The uso of the above method of obtaining closed expressions 

for tho distribution function is limited to very simple cases of 

linear lattices , whore the integral in (48) can be solved exactly, 

2 

For my reasonably realistic model tho dispersion relation, (jj (k), 
is canplicatod enough so that we have to resort to approximate 
methods. 


3, A pproxiaate Methods 

Tho most straightforward is the root-sampling method, in 
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which one generates a large number of froqtaencies by solving the 
secular equation at a largo number of uniformly distributed points 
in the first Brill ouin zone and then sorts these frequencies into 
sviitable ranges ( ifj , {jj ^ ^ so as to produce a 

histogram representation of g((iO ). The main objection to this 
nictViod is on the grounds of its being slow, since in order to 
nbtj.in a reasonably good sampling one needs a very large number 
of diagonalisations of the dynamical matrix. 

QQ 

Recently, Gilat and Dolling have introduced the so 
i;allcd "extrapolation method" which considerably accelerates the 

('.oitvrr.tional root sampling method. The idea behind this method 

A 

is to diagonalise the dynomical matrix relatively 

sriull number of evenly spaced points in the reciprocal space and 
then to find other solutions, for the eigenvffLues , in between by 
MC.^ns of a Taylor expansion about each such point for each eigen- 
value. 

Tho gradient of CO (q) needed for this extrapolation is 
crlculntcd numerically by imposing a small change ^ 
ith component of the wave vector. A now matrix & , representing 
tho chiingo in 'D(q) is constructed: 




(q 




/\0 

D 

jK 



( 50 ) 


Tho chmgo 



in the j-th eigenvalue is given by perturbation 


theory: 
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+ 


L 


I 


DJ < n a 

■ij‘s 


(51) 


^.t A 0 


where. - U U U matrix of oigernrectors 


A>0 


Of D . Hence we obtain 


( grart (/J ), 


g^i ^ 
g i*",* 


(52) 


aio dotorminos the value of CO . ((;^ + (S q) from the 

3 *** 

Tttlnr of CO (q) by \ising grad CO . as obtained above, i.e., 

j — 4 . J ■ -f 


^ i ( q + {^^ ) = CO (q) + grad OJ * A Q . 
J - - J - i 'j ' 


(53) 


in t.ho rictuvo]. calculation, the irreducible volume of the first 
'.r.'inoiun zone is divided into a uniform simple cubic mesh of 
!>''unt.i q . each noint being at tte centre of a small cube 
'.f.r ■':g’' 0 -it which extrapolation is carried out to a finite number 
of points. Bach cube has to bo properly weighted according to 
tk: nyi'inctiy of q^. The main point in this method is that 
.•o:-riuK-.hion of the gradient is much faster tlian re-diagonalisation, 
Gilat 'ind Dolling have n^ntionod that only tee first order term 
on t!io right hand sido of equation (51) need bo considered. At 
the aiim timo it should be noted that for q where has degenerate 
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rooos, tlie expression (51) is inapplicable and for such points 

r^O 

D h;xs to be rc-diagonalised to calculate C . ThG 7 have also 

t) 

pointc.d out that this method loads to relatively less accurate 
result;-, idierc tho surfaces of constant frequency display rapid 
v.’iriations. 

Qq 

More rc;ccntly Gilat and Raubenheincr have extended 
the extrapolation method for cubic crystals. This involves 
dividing the irreducible volume of the first Brillouin zone into 

uub.io "'.esh ,-and approxiraaiting the constant frequency surfaces 

* 

iur.id( every snail cuhj by a. sot of parallel planes perpendicular 
. Tho volume element dV falling between two 
uu h consecutive p3.ancs, describing OJ and + d uJ , 

in proportional to tho number of frequencies lying, between U> 

-iv! CO + d UJ , end for d lo 0 it can be approximated 

i,;,r 

dV = S dq , (54), 

V.*!’, n: 5 io tte cross-soctional area of a plane and dq is the. 
t f'dcknuay of tho volume clcnont. The distributxon function 
,;r given by 

g( uj ) « g (j, ^ J ) , (55) 

whi,n; tiiC' sumation is extended ever aXl eigenvalues and cubes. 

The function g(j, Q, ;,60 ) distritation of IWqusneios 

obtained irm U) Me) 5y extrapolating throu^out the cube 


rM ia given by 
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S( U q.JCj) dCJ = C W S(i^) dU-for 

% 

iO. grad ifj. j ^ 

/CO'. + C*>- trad Oj, I 

^ ^ ^ ml j Iq^ 

= 0 elsewhere, (56) 

where W is a weight factor associated xd.th the symmetry of q and . 

T-C — 

C is an arbitrary momber which is constant throughout the entire 
BrillcTuin zone.ij^j-is tlie distance of any plane from q and is 

the ma:cimum value of 05 .- inside a cube. The frequency {jj and 
the distance uS- are related by 

■ w = to, (q^) +' U> grad U). . (57) 

J 1 J £ = £0 

'I 

Gilat and Haubenheimer have given the expression for the cross- 
section area, in terms of the direction cosines of the gradient 
and the distance {jQ- , Ir. this method the entire Brillovdn zone 
and not just a finite number of selected points contributes to 
sioj) and this is precisely the advantage it has over the 
extrapolation method. Tlie errors, due to degeneracies and due 
to neglecting higher order contributions to the gradient, that 
occur in the extrapolation r-icthod persist. There are additional 
errors due to the assumption that constant frequency surfaces 
are parallel planes within the cube. 

/n other metliod which has had some popularity is the 
’'•'.oment method. TMs involves expansion of g((jkJ) in a series 
of Lefjendre polynomials in which the coefficients of the successive 
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terms are linear combinations of the nonients of the dynamical 
matrix, g ( OJ ) is written as 

00 

g( UJ ) = a p ( to/cJ ) , (58) 

n=0 ^ W t 

wliorc U>L is the maximim frequency and 

1 . 

gC IjJl X) P^(x) dx . ' (59) 

Since g( C*J ) is an even function of <LaJ , only the oven 
coefficients, are nonvanishing . Therefore we have 

.1 

(4n + 1) . f 

° ~T 3 — 'J 

L 


a 


(2n + 1) 


n 





2n 


= U, 


2n 


(60) 


where is related to the 2n-th moment (/^j^) • 

A 

iii— . = 

3 S 




to: 


’ “ vr-„ JZ 


(61) 


k 


The moment method has been applied by Montroll and 
collaborators to the evaluation of frequency spectra of a 
two dicionsional square lattice and of simple and body-centred 
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cubic lattices. 'The chief advantage of this nx;thod over the 
root sfjnpling methods is that one gets an analytic expression 
for g( Co ). One of the major deficiencies is tha,t it is 
inaccurate in the vicinity of a singularity due to slow convergence 
of tic Legendre polynomial expansion, 

93 

An improved method due independently to Lax and Lebowitz 
94 

t.nd to Rosonstock takes into account the presence of singularities. 
This method assumes tlmit it is possible to find the singular part, 

(x) , of the distribution function, g (x) is subtracted from 

S 

g(x) i;nd the remainder is approximated by a linear combination 
of n Lc;gGndre polyncainls . The first n-2 coefficients of the 
exparision are detcrm.incd by requiring tha.t the first n-2 moments 
of the approximate function agree with the exact ones. The 
remaining two coefficients are fixed by the value of the distri- 
bution function at the end points. This method is very efficient 
in the sense that only a few moments yield quite an accurate 
bchavio\ir of the distribution function. Despite the great amount 
of work on the location of critical points, in partic-ular by 
van Eovo"^'' , Rosonstock and Phillips , no general scheuB has 
boon found that will ensure that all the critical points have 
boon found. The applica.tion of this raethod, thorefore, is 
restricted to analytic models as one requires the knowledge of 
location raid shapes of singularities* 

98 

■ C, Isenbcrg has introduced a new method for the 
c;alculi.;tion of a large numbor of moments for crystals with short.^ 
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range interactions. He has proved a theoren that states: 

If /^2n (W) is the 2n~th jaoaent of the spcctrm for 
a finzto lattice , with cyclic boundai^^ conditions , in which the 
intci-actions extend to T neighbours, then 

= /^2n- ^ ^ ' provided n ^/T (62) 

3 

where H is the nmbc-r of atoms in the finite lattibe. The 
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use of the theorem aloi:ig with the relations 

S + t) — 0 
1 ‘1 


"j ^ i-i Vi " ”2 * " ' ° 


s 

n 


S + 
n-1 




+ 



for n ^ j 

(63) 


saves the labour of calculating moments considerably. Here 
is the nth trace of the jxj d,Yncmical matrix for a point in 
the first Brillouin zone and p.'s are the unknowns ttiat can be 
de to mined from the first j equations. The calculation 
consists in ovaluntli^ the first j traces for a point in the 
Bril3.ouin zone. The subsequent traces are obtained from equations 
(63). This is repeated for all the allovx;d wave-vectors. The 
corresponding values of S 's arc summed up, after siaitably 
weighting them according to the position of the point in the 
zone, to give us finally the first n even moments on the finite 



lattice and hence on the infinite lattice because of equation (62); 


Mention nay nlso be nade of Houston* s method which has 
onjo,'\'’C'd some popularity because of providing veiy good results at 
low frequencies, Ihis method uses expansion of the function 

distribution function per unit solid angle 
for tie Jth branchy in toms of Kubic harmonics. 


g 


OO f 

) K^( © , 9 ), 


(64) 


m=0 


where prime on the summation sign indicates the omission of the 

term crorre spending to m = 1. Now along the symraetiy directions 

the secular equation can be solved analytically for = *j 

© ,9) and for simple enough models this equation can be ■ . 

inverted to obtain k and hence g.(oj, 0 ^ ^ ) which is simply 
2 ^ 

(1/3- oV) k-^ (OJs Q f f ) dk^(uj , 0 ,(p )/dUl , The number of 
toms that we should retain in (64) is determined by the number of 
directions along which g can be obtained in an exact 

manner. The coefficients a (6J ) are evaluated by solving a set 


of simultaneous equations ; The total spectrum is then given by 

rr itt 

gj ( UJ ) = J sin 0 d e j gj (u), Q f f ) d(p . 

0 


(65) 


Houston's method has a serious drawback of introducing spurious 
sqtiare root singularities into the spectrum. Moreover, ihe derived 
frequency spectrum will converge to the actual one only if 
€( 0 J/&f^)ds evaluated along a very large number of directions, 
This itself indicates the ir^jractibility of the method. 
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Hwang has suggested interpolation between syrainetry 
lines or planes to obtain frequency contours within the Brillcuin 
zone. The rest of the procedure is aane as in Houston's nethod. 
The only calculation carried out by Hwang was for the tv;o 
dinensional Bowers-rRosenstock model. This method when combined 
vdth an exact handling of singularities gave excellent agreement 
v.dth the exact spectrum. 


Still another method that requires the use of contours of 
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constant frequency is due to Hwang and Wang , A few icroortant 
contours of ccsnstanf frequency on the surface of the BriUmin zone 
are calculated. The contours on the inner planes are traced by an 
interpolation s cheme . Urs volume enclosed by the surface of 
constant frequency is computed by numerical integration and fitted 
in a form 


2^/2 \ 3/2 2n 

■ X ' 

where the first tem is put according to result of singularity ana- 

95 96 2 2 
Ij’sis ^ . Finally G(CJ ) is obtained cn differentiating V(^ ) 

■ 2 ' 

vdth respect to OJ . This method was applied to the tw dimensional 
and three dimensional, simple cubic lattices and an excellent agreement 
with the exact frequency distribution was obtained. 


In short, these methods can bo classified into two main 
groups, tie sacpling methods and methods involving series eopansxons. 
The sampling. methods enable a histogram approximation of the fre- 


quency Spectrum, The advent of digital computers has led many 
workers to carry out e.xtensive; sampling calculations.. The .^ccur 7 
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IS greatly improved due to modifications by Gilat and Dolling, and 
Gilat and Ro-ubenheimer, These methods require the spectrum to be 
stored nmerically in the form of a histogram. Therefore numerical 
techniques have to be applied to obtain the thermodynamic properties 
iuid the accuracy of the final result is not easy to, assess. This 
makes it difficult to compare related models, itoong the series 
expansions, Houston’s modified methods, due to Hwang, and Hwang 
find Hang, have the disadvantage of requiring the knowledge of 
critical points and the constant frequency contours which are 
difficult to obtain even on the faces of the Brillouin zone, except 
for vciy simple models , The moment method, as such, is slow due to 
the poor convergence of the Legendre polynomials, and though the 
convergence can be improved in a manner suggested by Lax and 
Lebowitz, the computation problems associated with improving the 
convergence are formidable because of the need to know the location 
of critical points, Isenberg’s method permits us to calculate a 
largo number of moments on lattices where only short range inter- 
actions are present, A new series expajision method is sx^ggeeted 
in th,c following. 


Tlie Fourier Series Expansion Method 


The present method involves the use of sine functions as 
the b.csis functions fexr the expansion of the distribution function, 
Porhf.'.ps the only reason Legendre pdLyncanials have been used tradi- 
tinnally is because the coefficients of expiansion can be evaluated 
in terms of mcments. The use of other orthomormal sets, such as 
sine functions, does not offer this advantage, HovK3ver, this 
difficulty is a minor one in the age of high speed computers. 
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There are natxiral reasons for the choice of sine functions, 

I 

he distributicn function of a three dimensional lattice is zero 

at the end points and is Tjsually continuous in betvreen, its only 

pecularities being socie Icinks leading to discontinuities in the 

derivative . The Fcurier series of such a function is highly 
103 

convergent and therefore leads to a more accurate and convenient 
repiv^sentation, is in the case of the Legendre polynomial represent- 
ation^ tho Fourier series representation of the distribution 
function provides a means of obtaining a series expansion of any 
additive function of the normal mode frequencies of the lattice, 
such as thermodynamic properties. 

WiD shall consider here the distribution function G(X) of 
2 2 

the quantity X = U) /CJ in one of , the branches . The results 
oon be extended to other branches in the acoustic and optical bands 
in an obvious manner. The dispersion relation for the branch under 
consideration can. be written as 


X = X(k), (67) 

where X(k) gives the functional dependence of X on the wave vector 

k. 



so that 
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A 


n 


2 


1 

f 

G(I) ain (n 71 X) dX 


kJ 

= 2/v^ J ain n Tf X(k)j 
where is the voluine of the first Brillouin zone, 


( 69 ) 


The evaluation of the coefficients A thus anounts to 

n 

solving the integral in equation (69) losing the explicit form of the 
dispersion relation in equation (67). For some simple models the 
integral can be done analytically, but, however, it can be evalu- 
.'■'.ted numerically for all, models. Equation(69) when substituted 
in oquation(68) leads to the required representation for G(X). 


5’. Applications 


The method described above is applied to calculate G(X) 
for a model where the integral for te evaluated analytically, 

■fe consider here a m-dimensional simple cubic lattice with nearest 
neighbour central and non-central force constants being equal and 

104 

each atom having only one degree of freedom . The dispersion 
relation in this case can be written as 


X(k) = 1/fem (m - cos - C'OS ....... ^c.os Ic^j^) , (70) 

wlicre is the ith component of the m-dimensional vector k. 

m 

Substituting (70) in (69) and noting that 7* = (27? ) , we 


have. 


n 


, m 

2/(2 TT ) 



ssin 


nfC/Zvx (m - c.os - d’os 
1 cos IEjj^) dk^ dkg .... dk^ (71) 


O 



m • 
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After Gs^anding out the integrand it can easily be seen that the 
integral vanishes for even values of n. The nonvanishihg coeffi< 



The normalised form of G(x) is given by, 


Gei(X) = 


L. (- 1 )" I Jo ( n )1 m) 

n=n L 2m ■'J 


I- (-J-) 


7T(2n+l) 


A ai .t l.TT 
V 2m / 


For the one and two dimensional cases, d.e,, for m = 1 and 
2, the convergence of the series in equation(73) is bound to be poor 
due to the occurence of infinities in G(X),'*'Even with 30 teiros in. 
the scries, the approjdmate distribution function for the one 
dimensional case was found to oscillate about the cuive obtained by 
exact calculaticxi. The oscillations in G(X) are duo to strong infi- 
nites at the ends. 


In the two dimensional case table 7 depicts the approximate 

G(X) with 30 terms in the series, along with the exact one calculated 
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by Montroll , The agreement is surprisingly good, except at 
X = 0.5 and X *= 0 and 1, considering the fact that no attempt has 
been made to handle the logarithmic singxilarity at X = 0,5 and non- 
zero values of G(X) at X = 0 and 1 in the manner suggested by Lax 
ard Lobowitz. The present method gives a sharp peak in place of the 



locarithExic Infinity at X = .5. 


TABLE 7 


X 


(X) 


Exact 


Coj^uted 


0 

, .6366 

. 0 

.05 

.6706 

.7117 

.10 

. .7096 

.6988 

.15 

. .7549 

. .7608 

.20 

, .8086 

.8160 

.25 

.8740 

.8685 

.30 

.9560 

.9730 

.35 

1.0650 

1 .0537 

.40 

1.2224 

1.2445 

.45 

1.4977 

1 .4887 

.50 

oo 

2 .6689 


The height of this perJc is expected to increase vh-th increase in the 
nmber of tenns in the series. It nay be noted that equation(68) 
\<70uld always lead, to zero valijes of G-(X) at the end points, l-Jhile 
this is unrealistic for the two dimensional case_, it is possible to 
extrapolate to X = 0 and 1 from the neighbouring points to get values 
close to the true ones, It vdllbeshown later that one can get a very 
satisfactoiy representation for G(X) of a lamellar crystal by 
maUfying this method in the manner suggested by Lax and Lebowit^.^ 


In the three dimensional case approximate G(X), taking T 
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15,30 and 40 terns, along vdth the one obtained by taking 80 terms 
in the series (which can be considered to be the exact one) is depi- 
cted in Table 8, for getting an idea of the convergence of the ser- 
ies. .'.'igure 4 shows a plot of as obtained by Bowers and 

n ^ , , 104 

Rosens took and Montroll , The points are the results of the 

Fourier -expansion method with 30 terms in the series. -With 15 

terms, G(X) agrees with the exact one wi-fchin 1^ except near the 

end points end the critical point at X = 1/3 where the error is 

4-5jg. Increasing the number of terms to 30 and 40 decreases the 

error but the agreement at the end points does not improve much. 

This is due to the inability of the Fourier series to reproduce 

1/2 

X behaviour of G(X) when X — ^ 0. In the following we shall see 
how this difficulty can be overcome when low temperature studies 
are made. 

6 . A Modified Functional Expansion Method 

For the one and two dimensional lattices the Fourier series 
expansion method leads to incorrect resul-ts. In the two dimensional 
case, however, G(X) obtained by this method agrees closely with the 
exact one over most of the region, so -that in those cases where the 
finer details of G(X) are mimpor-tant, such as the calculation of 
the thermodynamic properties, G(X) computed by the above method can 
bo used. To obtain a more exact form of G(X) one can include cosine 
functions in the series expansion to give nonvanishing valxies at the 
end points, and a suitable number of singular functiens to reproduce 
the desired singu3.arities , . For the purpose of illustration, we have 
chosen a graphi-be-type laminar hexagonal lattice with nearest 
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TABLE 8 

Frequency distribution function of a three 
diiaensional simple cubic lattice. 



Exact 

Computed 


with 

X • 

0 (X) 


o 



0 

15 terms 

30 terms 

40 teims 

.0 

0 

0 

0 

0 

0.025 

0.175 

0*183 

0.171 

0*169 

0.050 

0.254 

0.248 

0.257 

0.252 

0.075 

0.326 

0.329 

0.328 

0.324 

0.100 

0.393 

0.392 

0.393 

0.392 

0.125 

0.460 

0.458 

0.458 

0.460 

0.150 

0.532 

,0.537 

0.530 

0.530 

0.175 

0.605 

0.598 

0.608 

0.605 

0.200 

0.689 

. 0.696 

0.691 

0,687 

0.225 

0.778 

0.774 

0.779 

0.778 

0.250 

0,887 

0.887 

O.'SSl 

0,884 

0.275 

1.011 

1.019 

i:oii 

1*013 

0.300 

1.183 

1.162 

1.188 

1-.181 

0.325 

1.446 

1,500 

1.479 

1.464 

0.350 

1.727 

1.745 

1.745 

1.731 

0.375 

1^727 

1.725 

1.731 

.1.723 

0.400 

1,720 

1.718 

1.720 

-1,718 

0.425 

1.717 

1.722 

1.715 

1.715 

0.450 


1.711 

1.714 

1.712 

0.475 

1.712 

1*717 

1.714 

1.711 

0.500 

1.713 

1.710 

1.714 

. 1.710 


1 
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* i Pi R 

neighbour interactions* The dispersion relations for this 
nodel are 


^ ± 1/6 4 cos 0 (cos <0 + 


cos 


<|j) + 1 J 




( 74 ) 


"here 0 and. ^ are proportional to the rectangialar con^onents of 
the wave vector in reciprocal space and range from 0 to 
•~,:)sitj.vG and negative signs correspond to the two branches. GCX) 
has boon evalut?<.ted in a closed form by several authors (see e.g,' 
ref. 105). It is 


Cr(X) = — K(A/B) for 0 ^ X ^ l/Z 

yf* B N> 


Tt2 a 

and G(X} = G(1 - X) 


= JJkjzML K (b/a) for 1/3 /C' X / 1/2 (75) 

for l/Z 

for X <^0 cr X^ 1, 


= 0 


where 


.cl 


3 2 

3X (2 - 3X) , B = 3 (1 - 2X) 


and 


Tr/2 2 2 -J-/2 

K(Z) = J (1 _ Z sin ^ ) d ^ . 

0 

G(X) contains logarithmic singularities at X = l/3 and X = 2/3, 


To obtain the distribution function by the modified Fourier 
expansion method, to have included two extra terms, B^ cos .JT f and 

In I f - 2/3 j , the former to give a nonzero value at the 

^ » * 

origin and the latter to reproduce the singularity at the point 
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r i 2 2 

? =2/3, whore f is LjJ "the acoustic branch and is 2X. 

T,!ic d:iotiubution function for Jhis branch is 

G(f) = cosrrf + In I f _ 2/3 I +X sin n n f (76) 

n=l 

In .ua, 7 ' be mentioned that a pidor knowledge of B and B is not 
n',.'eded.. B^ and are calcula.ted by solving p + 2 sinultaneous 

linear equations. It has boon found, by trying siiccessively 
incroas'.ing values of p, that only seven sine functions are suffi~ 
cient to yield a satis factoiy agreement with the correct G(x) given 
by eq\i£,tions (75). The computed values of the coefficients are 

B = .302 j B = ~ .275 * = 8.2753d0 ; 

1 2 ^ 

•-2 —2 -2 
Ag = ^.909x10 ) ■ Aj = - 1.6803d0 ; = -1.928x10 ; 

-2 “3 

Ag = 1.962x10 ; Ag = -8,772x10 ^ kj = -6,530x10 . 


The distribution function for the other branch can be obtained 
by reflexion of G(f ) given by equaticai(76) about f = 1,0 axis for 
t is model. In figure 5, tlie solid curve is the exact result from 
e.'Uiations (75) and the points represent the values obtained by this 
nodifiud mothal. 


In exactly the 'same manner, in the three dimensional models 
a bettor agreument between the approximate G(x), and the exact one, 
at tlic end points, can be achieved. The behavioiir of G(X) at the 
end points is known to be^*^^ 


g(x) oC X 


1/2 


f or X -> 0 

1/2 


(77) 


<?C (1 - X) 


f or X 1 
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X 


FIG. 5 G(X) FOR A GRAPHITE - TYPE TWO DIMENSIONAL LATTICE 
WITH NEAREST NEIGHBOUR CENTRAL FORCE INTERACTION. 
THE SOLID CURVE REPRESENTS THE EXACT CALCULATION 
DUE TO ROSENSTOCK, AND THE CIRCLES SHOW THE RESULTS 
OBTAINED BY THE PRESENT METHOD. 
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1/2 l/2 

how if wo introduce a term X (1 - X) in the series expansion^ 
G(X) will still vanish at the end points but will behave in the 
neighbourhood as required (77), 

- The Distribution Function for some F.C.C. and B.C.C, Lattice 

bb proceed to apply the unmodified Fourier scries expansion 
moohod bo compute the distribution functions of some f.c.c, and 
b.c.c. lattice models, the modification being not necessary as shall 
br. pointed out in the next section. For f,c,c. lattices we assume 
a central force model. For b,c,c, lattices wo take a model in which 
each atom interacts with its nearest neighbour through central and 
angular forces and with its next nearest neighbours throu^ central 
forces c For both typos of lattices the interactions have been 
liiaited upto second neighbours. The secular equations can be 
written as follows ; 

( a) F X. , C- , lattice 

1 2 2 
2 — tp(C2 + Cj) + 2 flT" S-j^ ^ ^1^2 ^1^3 

2 2 

2 - ^ *” ^ ^2^3 

2 v2 

S„S, 2- C (C. + Co) + 2o-S_ - A 

x3 #0(5 o JL ^ V 


= 0 . (78 } 
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(t>) B.C.C. Lattice 


j(l+2 f 

J ^ ^ \2 

i 

I 



1(1- 





12 3 


(1-P)V2C3 (1- P )S3_C^S3 


(1+2 p ) (1-03^0203) (1- P ) 32 S 3 C^ 

+ 1 <r' ^ 



(1- {3 )S2S^C^ 


(1+2 p )Cl-CiC2C3) 

+ I tr- s^- >? 


(79) 


whore C = cos k - S. = sin k. , k. = q. q ,2a = lattice constant 
j 3 J J J 0 

2 

M CjO /2c<. , for F.C.C. lattice 


3M U> /80C , for B.C.C. lattice. 


and )\ = 


Horo M is the mass of the aton. in i±ie lattice, and (j~~ are 

iiho nearest neighbour noncenti^ and next-nearest neighbour 
central force constants expressed in rants of the nearest 
noighbotir central force constant . 1, J'th rectangular 

coiinpcnont of the wave vector q, 

V 2 

Solving the cubic equation for A ^ dispersion 

relation for each frequency branch. For each branch g(x) can 
b ) evaluated by using equations (68) and (69) , taking care to 
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ensure that the integration in equation (69) is done through 

the Brillouin ssone corresponding to the lattice model. In these 

cubic lattices , the symmetry of the constant frequency surfaces 

can be exploited by doing the integration through only 1/48 of the 

first Brillouin zone, which in the F.C.C. lattice is the part of 

the zone enclosed in the trihedral angle defined by (001), (101) 

and (111) directions, and in the B.C.C. case is the part enclosed 

x^ithin the tetrahedron bounded by the planes k = 0, k +k = rr, 

2 1 o . 

k^ = kp and k^ = k^, with their ranges defined by 


and 0 <“3 4.Tr . 

The integration required in equation (69) was carried out by a 

CDC-3600 electronic computer.. For the F.C.C. lattice we ha-ro 
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used the procedure given by Overton and Dent to evaluate the 
frequsiicios in a particular branch in the region of integration. 
For the B.C.C. lattice the volume integral is expressed as a sum 
of two three dimensional integrals: 


f 


n/2 k k. 




r 

r" r"- r 1 

sin 

rnTT,)tCk)] d\= ' 


dk^ sin n TT X (k) 


i 


0 -^0 
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Both the iiitegrals on the right hand side of (80) are written 
in the same program, with the machine instructed to change the 
upper limit of the variable from to as soon as 

the value of k^ exceeds TTA. For the purpose of integration 
the whole volume of the irreducible part of the zone is divided 
into cxibes of edge length TTA^O along each axis. Thus by increa- 
sing kj, and in steps of TT/lOO up to their respective 
niaximum limits, the machine calculates the value of the integrand 
at tho centre of each'' cube and adds then by storing them at the 

same location. The final result is obtained by mtiltiplying it 

g 

with the volume of the cube ( tT /lOO) . 

We have calcxiLated the first fifteen coefficients of 

tho function sin (ntT^) in equation (68) for eleven values ofo^ 

ranging from -.25 to +.25 in steps of ,05, for the face-centred 

cubic lattices - and five values of ^ , frem - cr*"/5 to + cT'/S 

in steps of <r“/l0, corresponding to each of the eleven values 

of cr* , from 0 to 1 in steps of .1, for the body-centred cubic 

lattices. Table 9 contains these coefficients for three values 

of if for tho face-centred lattices. Tables. io(h) ijnd 10(b) contain 

tho coefficients for three values of |3 'with = 0 and three 

values of (f* with p =0, respectively for the body-centred 

cubic lattices. Here L, T^ and T^ refer to the longitudinal and 

two transverse branches respectively. The dimensionlass variable 

X is the ratio , so that the range of X is fren zero 

^max 

to one for each branch. The coefficients have been so noimalisod 
that the area enclosed by G(X) for each branch is one, when G(X) 
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Face - centi«d cubic lattice ^ 
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cr 

= -0.1 


<r' = 0 


• 

o 

1! 


n 

L 

Tl 

T2 

L 

7l 

% 

L 

^1 

% 

1 

1.146 

1.675 

1.473 

1,051 

1.723 

1.558 

1.123 

1.734 

1.612 

2 

-1 .172 

0.913 

1.118 

-1.142 

0.704 

1.110 

-1 .256 

0.682 

1.168 

3 

1 „112 

-0.407 

0.234 

1.147 

-0.612 

0.094 

1.217 

-0.680 

-0.052 

4 

-0.886 

-0.262 

0.138 

-0.982 

-0.161 

-0.134 

-0.934 

-0.177 

-0.440 

5 

0.574 

0.211 

0.106 

0.751 

0.314 

-0,199 

0.551 

0.381 

-0.298 

G 

-0.229 

0.005 

-0.170 

1 

O 

-0,060 

-0.197 

-0.118 

-0.035 

-0.012 

7 

-O.C'GS 

-0.142 

-0.118 

0.156 

-0.146 

0.102 

-0.227 

-0.170 

0.295 

8 

0.309 

0'.002 

0.046 

• 0.104 

0.091 

0,143 

0.'443 

0.120 

0.168 

9 

-0.435 

0.015 

0.008 

-0.276 

0.075 

-0.0C2 

-0.478 

0.095 

-0.111 

10 

0.460 

0.017 

0.089 

0."364 

0.025 

0.128 

0.373 

-0.029 

0.025 

11 

-0.378 

0.097 

0.173 ■ 

-0.'350 

0.024 

0.098 

-0.168 

-0.011 

0.087 

12 

0.236 

0.037 . 

-0.029 

0.266 

-0.053 

-0.173 

-0.046 

—0.040 

-0.185 

13 

-0.057 - 

■0.034 - 

-0.099 . 

-0.126 . 

-0,023 

-0.071 

0.216 

-0.007 

-0.102 

14 

-0.103 

0.016 

0.120 . 

o 

-1 

0.046 

0.144 

-0.289 

0..081 

0.229 

15 

0.2.'". 5 

0.013 

0.110 

0.157 . 

-0.-016 . 

-0.016 

0.269 

0.005 

0.104 


I .refers to the longitudinal branch, and and to the 


transverse branches. 
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TABIE 10(a) 

Body-centred cubic lattice , =■ Q 




*= -0.05 j 

^-0 

f 

B= 0.05 

n 

L 


To L 

T. T„ 

L 


T.. 

1 

2 

1 2 


1 

2 

1 

1.081 

1.072 

1.703 1.066 

1.075 1.703 

1.056 

1.082 

1.702 

2 

-0.994 

0.940 

0.153 -0.990 

0.945 0.150 

-0.988 

0.952 

0.147 

3 

0.925 

0.839 

-0.614 0.921 

0.840 -0.616 

0.919 

0,841 

-0.617 

4 

-0.718 

0.619 

-0.070 -0.717 

0.616 -0.070 

-0.714 

0.611 

-0.068 

5 

0.546 

0.502 

0.519 0.544 

0.495 0.515 

0.541 

0.486 

0.519 

6 

-0.361 

0.401 

0.079 -0.365 

0.393 0.076 

.,.0.365 

0.382 

0.074 

7 

0.242 

0.368 

-0.347 0.253 

0.358 -0.349 

0.260 

0.346 

-0.350 

8 

-0.141 

0.306 

-0,070 -0.164 

0.296 -0.070 

-0,181 

0.283 

-0.068 

9 

0.100 

0.265 

0.317 0.134 

0.253 0.318 

0.161 

0.240 

0.319 

10 

-0.079 

0.207 

0.070 -0.124 

0.197 0.068 

-0.160 

0.184 

0.067 

11 

0.095 

0.183 

-0.234 0.146 

0.175 -0.236 

0.186 

0.164 

-0.238 

12 

-0.107 

0.155 

-0.058 -0.161 

0.150 -0.059 

-0.201 

0.142 

-0.058 

13 

0.127 

O.lStL 

0.224 0.177 

0.148 0.224 

0.212 

0.141 

0 .224 

14 

-0.124' 

0.136 

0.053 -0.168 

0.134 0.052 

-0,197 

0.129 

0.052 

15 

0.115 

0.133 

-0.178 0.152 

0.134 -0.180 

0.173 

0.129 

-0.180 


L refers to the longitudinal branch, and and Tg to the 


transverse branctes , 
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TiBIE 10(b) 

Body-centred cubic lattice, ^ = 0 



= o,z 

! 

<r 

- = 0.4 



/T' = O.E 

5 

n 

L 



L 

.^1 

"2 

L 

^1 

T 

2 

1 

1.065 

1.530 

1.630 

1.103 

1.629 

1.591 

1 

.131 

1.686 

1.531 

2 

-1.145 

0.948 

-0.415 

-1.202 

0.712 

-0.524 

-1 

.242 

0.397 

-0.589 

3 

1.133 

0.223 

-0.228 

1.167 

-0.175 

-0.036 

1 

.176 

-0.476 

0.129 

4 

-0.934 

-0.244 

0.329 

-0.912 

-0.426 

0.20s 

-0 

.864 

-0.330 

0.025 

5 

0,688 

-0.256 

-0.023 

0.593 

-0.133 

-0.076 

0. 

,489 

0.178 

0.034 

6 

-0.398 

-0.154 

-0.088 

-0.243 

0.120 

0.089 

-0. 

.109 

0.271 

0.088 

7 

0.159 

0.026 

0.110 

-0.024 

0.224 

-0.059 

-0, 

.144 

0.054 

-0.166 

8 

0.031 

0.121 

-0.080 

0.201 

0.083 

-0.038 

0! 

,271 

-0.178 

0.132 

9 

-0.135 

0.129 

0.064 

-0.254 

-0.067 

0.169 

-0. 

.261 

-0.097 

0.012 

10 

0.180 

0,036 

0.052 

0.230 . 

-0.110 

-0.122 

0. 

,189 

0.091 

—0 ,054 

11 

-0.161 

-0.020 

-0.107 

-0*147 . 

-0.005 

0.013 . 

-0. 

,090 

0.127 

0.056 

12 

0.117 

-0.034 

0.048 

0.061 

0.061 

0.058 

0. 

.017 

-0.021 

-0.056 

13 

-0.056 

0,002 

0.002 

0.011 

0.053 

-0.038 

0. 

.031 

-0.081 

0.088 

14 

0.012 

0.012 

-0.096 

-0u045 . 

-0.007 

0.038 

-0. 

.052 

-0.018 

-0.053 

15 

0.014 

0.012 

0.029 

0.053 

-0.020 

-0.02-8 

0. 

,070 

0.063 

0,000 


L refers to the longitudinal branch, and and Tg to the 


transverse branches 
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A 


FIG. 6 FREQUENCY DISTRIBUTION FUNCTION G(x) OF A F. C. C. 
LATTICE MODEL WITH a= -0. 10. THE SOLID CURVE 
REPRESENTS THE PRESENT CALCULATION, AND THE 
CIRCLES SHOW THE RESULTS OBTAINED BY OVERTON 
AND DENT. HERE = mJ- / 2a. 
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FIG. 7 FREQUENCY DISTRIBUTION FUNCTION G(f) OF VANAU.IUM 
TAKING THE NONCENTRAL FORCE MODEL, f IS THE 
NORMALISED FREQUENCY, co/cOj^ax' 
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is expressed in a Fourier series with the number of terns in the 
series limited to the maximum value of h given in the tables* It 
may be mentioned that the machine time needed for computing the 
coefficients for one model was about one minute* 

The frequency distribution function G(^) for 57 ^= - .10 
in a face-centred cubic lattice is compared witti the calculations 
of Overton and Dent in figure 6, . The agreement, in general, is 
quite satis factor^'-, . 

The distribution function of vanadiuum has also been 

107 

computed by using the elastic constants given by iSlers and is 

plotted in figure 7, Oon results are in qualitative agreement with 

108 

those obtained by Hendricks et al, iising the root sampling method. 
® • Be mark on the Fouiier Expansion Method 

The method for calculating the frequency distribution 
function of a crystal presented above has two distinct advantages 
over the conventional root sampling method - a very considerable 
economy of time and numerical work, .and yielding simple analytic 
expressions for the distribution function which can conveniently 
be used in computation of quantities that require the distribution 
function. The accuracy of the method can be enhanced by increasing 
the number of terms in the series and at the same tine using a 
finer raesh for the integration invcO-ved in equation(69) . This 
method does not reqmre a close examination of the topology of the 
constant frequency surfaces particularly the points of degeneracy 
where two or more branches cross* “Rie methods of GLlat and Dolling 
and Gilat andHtnxbehheimer reqiaire prior knowledge of these points 



..84 


of dej^eiJoracy in k space. 

!Ilie convergence of the Fourier series is mch better than 
the corresponding series in Legendre polynomials. However the 
ntanber of terms one must retain to get a satisfacto:yreprcsentation 
of G(X) in a given case is best determined by trial. The extreme 
complexity of tile dispersion relations makes it difficult to 
establish rules on this point. Wo have found thft 15 to 20 terms 
give quite satisfactoiy results. The convergence of the Fourier 
series for the longitudinal branch is poorer than that for a 
transverse branch because of the peaking that occurs in the former 
in the high frequency end of the spectrum. 

Since a finite Fourier series can never reproduce the 
critical points exactly, the question of how closely can the 
location of the critical points be determined by this method is 

of some interest. The critical points in G(X) are usually such 

// 

that they lead to ^ -function singularities in G (X) , If one first 
ovaluatosG(X) with sufficient accuracy i.e. by keeping an adequate 
number of terms in the Fourier series such that convergence has 
become obvious in the form -of the smallness of the last few coe- 
fficients, G^(x) can then be evaluated in the form of the series, 

gI^X) ■ = -7T £I ^ sin(nTt X). 

■DvB to the poor convergence of the Fourier series at the 
critical points, as we have seen in the case of a three-dinensional 
simple cubic model, we shall generally get a peak of finite width^ 
instead of g -function type singularity at such a point, Thea?efore, 
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the nost one oan say is that the critical point lies within the 
width of the prorainant peak* 

It renains to say something as to the law of decrease of 
the successive terras in the Fourier series eapansion. It is evident 
at once that the values of the coefficients An must u].tinately 
vanish indefinitely as n increases, owing to the more and more rapid 
fluctuation in the sign of sin(n ty X) and the consequent more 
complete 'cancellation of the various elements in the definite 
integrals for the coefficients. Generally if a function f(x) is 
continuous and its derivative has a finite nun'ber of isolated 

103 

discontinuities {in a period), the convergence is ultimately as 


1, -JL_ 

22 






The distribution function for any physical three dimensional lattice 
is known to be continuous, A finite number of discontinuities occur 
only in the derivative of G(X), Therefcre, the cornu rgence of the 
sine series expansion for G(X) is the same as stated above, 

121 

Recently Mahanty has suggested a method to evaluate the 
lattice Green functions using a Fourier series expansion of their 
imaginary part which is related to tho real part by a simple disper- 
sion relation. The method has proved to 'be much faster than the 
direct integration of the Green function integrals and has the 
added advantage of yielding analytic expressions for the real and 
imaginaxy parts of tho Green functions. 


Having established the usefulness of tho Fourier series 
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expansion method, we proceed to apply it to the case of diamond in 
the following chapter. 



CHAPTER IV 


APPLICATION OP THE FOURIER SERIES EXPANSION METHOD TO DIAMOND 

1 , Frequency Distr ibution Function 

The elerients of the dyneanical natrix for a Born-von Kanaan 
nodea with interactions upto second neighbours were obtained in 
Chapter II, The values of the six force cons-tants i£s obtained by 

73 

different methods have been shown in Tab3e 6. Smith was the 

first to calcifLate the frequency distribution function of dianond 

(two neighbour Born-von Karman model) with the force constants 

which she had determined from Raman frequency and the elastic 

constants and C^^. The frequencies of the normal modes in the 

(100) and (HI) directions with Smith’s force constants have been 

found to be novhere near the experimental values from neutron scatt- 
12 

ering results , Though the force constants calculated by Saint 
J,anes do not yield better agreement, the distribution function has 
been obtained for each set of force constants calcxiLated by him 
and the force constants calculated in Chapter II, in order to see the 
effect of variation in force constants. 

The integration in equation(69) is done through the irr- 
educible part of the first Brillouin zone which is the sane as for 
a f.c.c. lattice and hence the sane procedure of integration is 
rdloptod. The dynamical matrix now is of order 6 and has complex 
elements and, therefore, the diagcnalisation consumns more time. 

It may be pointed out that one can do away the need to know the 
naxinm frequencies in the various branches , vhich probably appears 
ossenti;:! from the analysis of the Fourier series expansion method. 

R7 . 
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This is dito + 4. 

0 the fact that one can always consider the distribu- 

io 

Action, with zero value ^ to exist outside tlie eillowed band. 

Therefore •i-p ^ > • 

* uj IS normalised with respect to any arbitrary frequ- 
Gncy + 1 • 

uaking care that is not less than the maximum 

quency of the allowed baoid, then the Fourier expansion nethod 
will load +A , 

o a function which represents G(X) inside the allowed 

C'Xlci Vi • 

oscillations, with small amplitudes abo^it the X-axis 
"de the band. g(x) can be truncated at a point where these 

Oscill.,-tinr,o ^4. 4. 

^ --tart. This fact was made use of in the case of 

diamond .• 


'^he initial guess for O-J is the frequency of the optical 

photlCiJi aV 14. 

^1=0, which is 2.51 x 10 rad/sec. The arbitrary 

q'ucncyj therefore, is chosen to be 3x10^"^ rad/sec. to be on the 
3-afer aide. 


coefficiGnts for the first 30 terms in the Fourier 
given in Table 11, and the corresponding distribution 
^ are plotted in Figure 8, 


chape of the distribution curve, in each one of the four 
> is .very much different from the one obtained by 

Ih.iith 

i '-'O much so that we get four peaks in the spectrum cs 
-v)-'-innt only two pools occuring in the later case. Smith had consi- 
■■'d only central interactions with the second neighbours and calcu- 
™ spectrum by root sampling method by talcing only 29 points 
irreducible part of the zone. It nay not be irrelevant to 

ntiin that the spectrum for Ge and Si, which was calculated by 

Colo {ind Vt 11 0 

tvinQkB , also shovrs four peaks. 
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TABLE 11 

Coefficients of the 30 terns of sine series 
(force' constants used are frai table 6) 


11 

I 

II 

III 

■ IV 

1 

4* 

25il68 

26.045 

25.101 

24.681 

o 

1 . 

5.897 

3.724 

7.705 

4.771 

7 , 

-1.912 

-3.641 

-2.601 

-2.590 

4 

5.916 

6.097 

4.820 

12.006 

5 

-0.035 

-1.156 

2i463 

3.724 

6 

-1.428 

-2.756 

0.273 

-5.808 

7 

,1.142 

1.250 

0,210 

3.326 

8 

-0.171 

0.100 

-1.077 

2.536 

9 

0.491 

-0.311 

0.892 

-5.225 

10 

-2.904 

-1.-397 

-1.450 

-0.904 

11 

0.950 

2.064 

-1.907 

3.573 

12 

5.007 

2 . 928 

3.277 

-0.256 

13 

»^1 . 845 

-3.092 

1.228 

-2.776 

14 

-2.417 

-0.248 

-1.970 

0.048 

15 

, 2,264 

2.651 

-0*222 

1,871 

16 

17 

0.985 

0.426 

-1.305 

-0,517 

0.660 

2.130 

-0.093 

-0.079 

18 

-0.908 

0.186 

0.860 

1,055 

19 

-0.997 

-0.100 

-2.161 

-0.723 

20 

1,639 

0.823 

0*146 

0.478 

21 

-0'^939 

-0.928 

lil37 

2.644 

22 

23 

2-.: 

-0.885 

2.052 

-0.705 

0.626 

1.511 

-1.759 

-0.467 

0.239 

-0.306 

-1.117 

-2.077 

1.760 

25 

-0,751 

-0.050 

-0.367 

•0.555 

26 

1.117 

1.298 

. 0.613 

-1.315 

27 

-0.048 

-0.596 

-0.542 

0.735 

28 

0,478 

-0,039 

0.265 

0.114 

29 

30 

-0.266 

-0.490 

1.209 

-1.400 

-0,520 

-0.085 

-0.600 

0 . 449 
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One inportan,t difference among the varioi:is cui^s in fig. 8 

14 

the valioe of the maximum frequency. For curve IV it is 2,55x10 
r.id/soc as againat the value 2,60x10^^ rad/sec for the other three 
curves* This is probably due to the finite value of ^ i^ 
fonnor case. It appears from Fig. 8 that the distribution ourve is 
very sonsitivo to the parameter ^ . 

2 , Specific Heat 

If Cj(X) is the distribution of the squared frequencies of 
normal modes of a crystal, the specific heat at constant volume 
mole is 




SNKgJ G(X) E ^ 

c 


ti (jj 

Kg! 


dX 


(81) 


whore X = , N = Avogadro»s number. Kg = Boltzmann constant, 

s- . _14 


h = Planck * 3 constant, T = Absolute temperature, ® 

rad/sc c and E(y) is given by 

, X 2 y . V ^2 

F(y) = y e / (e - 1) 


(82) 


G(X) is normalised so that the integral over the entire spectrum xs 

’xndty. The integral ih equation(8l) is solved numerically 

VAri'-rus values of T, using all the four distribution functions. The 

vi -rro been obtained 

cxT.'Orimcntal values of the specific heat of diamond na 

111 1^3 . .V- temperature 

by Nonint , Robertson et al. and Pitzer m the no 

r•J^gG 70 - 300 K, Comparison of the theoretical and expc 

. otrident that I and 

voaues of specific heat is made in figufe 9. It is evx 

1 + q yh6r^ ^-^ 

II sots of force constants do not give accuratd resuiPs, 

V + Tn close agreement 

III ,ind 17 sets produce values of the specific hear mx 



Cv 

(Cals /OK) 



T (0 K) 


FIG. y SPECIFIC HEAT OF DIAMOND USING THE VARIOUS PHONON 
SPECTRA FROM FIG. 8. 
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with the experiinetttal values . Figure 10 shows the relation between 
the experiiaental and- theoretical values (for model rv) of the 
characteristic temperature 

3, Localised Modes 

In pure diamond the first order dipole mcnont associated with 
any phonon is zero by symmetry-. The introduction ^ impurity 
with non~zero effective charge Isa-^s to absorption by creation of 
single phchoha. The presence of the impurity modifi®® the normal, 
modes of the perfect lattice, particularly in its vicinity. Hence 
the vibrational amplitude and therefore the dipol® moment associated 
with the impurity atom are not the same as those calciilated by using 
the unperturbed displacements. The absorption spectrum is thus 
characteristic of the impurity concerned. In the case of local 
modes, where the atomic displacements are confined to the vicinity of 
the impurity, this effect predominates in determining the absorption 
spoctrum* Dawber and Elliotb have treated ihe case of a 

substitutional charged impurity atom of mass M without any change 
of force constants. They show that fpr the case of monatomic cubic 
and diamond structure materials containing a mass defect, the fre- 
quencies of the perturbed modes are given by 

ih 

1 *=^Z f' G(X) dX (83) 

Jq Z - X 

wtere Z = CO / Cja gives the squared frequency perturbed 

nodes and X *= gives the squared ^ perfect 

lattice modes whose normalised density of state i^ G(X) ^'measures 
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c?iangc in mass, (M *- M)/M » Tlie solutions of equation(83) are of 
two kinds according as Z lies outside or inside the regicsn of 
unperturbed frequencies which range from 0 to a maxinura 

tl * 

value OJ . 

L 

For a light subs ti tut i'Xial impurity one triply degenerate 
mode nay lie outside the region of frequencies allowed in the perfect 
lattice. Such a node is called a localised node. We shall concern 
ourselves only with such nodes. The value of ^ above which a 
localised mode nay occur, called the ciitical mass, depends on the 
shape of the distribution functicsn used, 
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ingress et al. have analysed the absorption spectrum of 
Si - containing boron in great detail and have found good agreement 
with Dawber and EHiotJfs theory. 


We calculate here the localised mode frequencies for various 
values of the mass defect on the basis of the theory of Dawber and 
Elliott, using different frequency spectra of diamond. 


If we substitute the Fourier series for G(X) in equatian(83), 
we obtain 


ez r 




sinfnTTX') 


dX 


Z - X 


j^cos(nnZ)|^-Si(n7TZ) + Si(ntT(Z - 
- sin(nfr Z) ^Ci(n7XZ) + Ci (nj[(Z - , (84) 



.,96 


whore S±(y) and Ci(y) are defined as 

7 ■ ■ 

Sify) = f dt 


i(y) = r ■ 


Cl(y) = 


oo 

J 


cos t 


di 


(85) 


.^‘••■,''^' 13 ; tivc tables for the Si and Ci functiens exist (■ .ft, soc ref.?0}, 


The amplittodo of the defect atoEi,*^(0), is 
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X(o) I 


M 


f z 


^ ■G(x)dx _r 

I (z - if ^ 




( 86 ) 


For the Fourier seiles representation of G(X), equation(86) reduces 
to 

r . _ _ r 

lTTXL 


H 


I %(t)r 


,2 2 

^ Z 


£ i /ajanrl 

^ ^ • l Z - Xt 


+ n ff sin(n1T Z) 


X ^ Si (n 7T (Z X^) ) - Si (n TT Z)J 
+ n7Tcos(n7T Z) ^ Ci(n7T (Z ^ X^)) - Ci(n7TZ)J'-C 


-i 


(87) 


Figufe 11 shows plots 



and M 


^|9((0) I against . . 


It is at once clear that the curves (a) for the localised node frequeno- 
ier oie quite insensitive to the finer details of the frequent^ 
spoctnm. On the other hand the critical values of the mass defect 
paraiaetor,. €.q, are very nuch different for the different distribution 
curves. Though this statement is generally true, two different 
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distributioBE ctirves may lead to the sane value of ^ . viii-di in fact 

c 

is the case for curves I and TV, The values of ^ ■ for the distriba- 

C' 

tion curves obtained by xasing and IV sets of force constants 

are approximately ,23, .14, ,35 and ,23 respectively, Kt may be 

117 ^ 

mentioned that Angress et al. foxmd the value of tg as .25 
for silicon, using a partictilar form of phonon spectrum. 

The other set of curves, b, in figiare 11 represents the 
relative amplitude of the defect atom. These curves are con^ara- 
tively more sensitive to the form of the phonon spectrum for small ^ « 
The properties of the local mode become increasingly less sensitive 
to the phonon spectrum employed as ^ increases (which means as the 
impurity mass decreases). This is quite expected because the local 
mode frequency increases rapidly as £ increases towards unity and 
the vibrations be core more concentrated roimd the defect, and con- 
sequently the impurity atom vibratos in viidxially static surroundings . 

The absorption spectrum in diamond has been obsesrved by 

112 119 120 

Robertson et al, , Sutherland et al, and Smith and Taylor , 
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and the different spectra are compared by Elliott who has shown 
that only type I diamonds give a local mode peah at about ,17 eV 
2,6x 10 rad/sec) just above the maximum frequency of the allowed 
band. This cori^sponds to a mass defect just above the critical 

jf 

mass defect parameter which is ,23 for the phonon spectra I and IV 

9 

in figure 9-, This reflects the presence of Be isotope.. The 
ccnclusion must be taJcen as very tentative as it is btised on Bawber 
and Elliotts* formul.a(83) which is valid only for an isolated 
charged mass defect idthout any change in spring constants. 
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CONCLUSION 

In the thesis firesented above, firstly a new method based csn 
a Fourier series expansion Has been used for the evaluation of the 
freqtiency distribution function, it has been apj^lied to a two neighbour 
Bom-yon Karman model for faceucdhtred (with central force interactions 
only) and body-centred (with non-central force interactions with the 
iimt neighbours and central forde interactions wiidi the second neighbours) 
dubic lattices* It is found ihat 16 terms in the Fourier series are 
sufficient to obtain a reasonable lepresentation of the frequency distri<^ 
bution function. The first fifteen coefficients have been evaluated for 
these lattices with various force constant ratios. 

Secondly the lattice dynamics of diamond has been studied in 
detail. The simplified two neighbour Bom— von Karman model ( with the 
second neighbotiT' interaction being limited to only central forces) that 
was previously studied by Smilh is found to be in disagreement with the 
neutron scattering data obtaitied recently by Warren et. al* For this 
reason the angular and torsional forces haVe also teen included in the 
interactions with ihe second neighbodre* ^ch a model contains only five 
parameters. The three elastic constants and the Raman frequency were pUt 
as constraints for the evaluation of the force constants.- Only one para-: 
meter ^ was allowed to vary to get the best fit with the measured fre- 
quencies at the zone boundaries in the (100) and (111) directions. When 

76 

the elastic data of Bhagavantam and Bhimasenachar is used,.- the value 
= ,044 is found to give a reasonably good- fit with the measuied 
dispersion curves.. The deviation is generally within the experimental 
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